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CHANNEL-SLOPE FACTOR IN FLOOD-FREQUENCY ANALYSIS 


Manuel A. Benson,! A.M. ASCE 


ABSTRACT 


Annual flood peaks in New England have been related to many hydrologic 
tors. The main-channel slope has been found next in importance to 
1inage-area size. The slope for that part of the main channel between 85 
110 percent of the total distance above the gaging point provides the best 
rrelation with flood magnitudes. 


Many agencies and individuals are studying the causes of floods and their 
ation to known meteorologic and physiographic conditions, and are develop- 
generalized flood-frequency relations for the purposes of hydraulic de- 
n, flood insurance, flood-plain occupancy, and other uses. 

Data for New England have been used in the first phase of such a study by 
U. S. Geological Survey. The records of runoff and precipitation for this 
ion are aS numerous and as long as records for any other region in the 
ited States; topographic quadrangle maps covering the entire area are 
lable; historical flood data are available, dating back to the earliest 

es of settlement. 

Peak discharge records of 170 gaging stations were used in the analysis. 
se records represented all the stations in New England having at least 10 
rs of record, on which the effect of regulation on annual peak discharges 
3 known, or judged, to be negligible. These stations were located through- 
the region, though not as densely in the northern part. The range in the 
= of drainage areas was from 1.64 to 9,661 square miles. Of a total of 
Station, 3 stations drained areas less than 10 square miles; 31 stations 
ined areas less than 50 square miles; 65 stations drained areas less than 
Square miles; and 106 stations drained areas less than 200 square miles. 
n practice, there are many ways of converting a set of flood data at a 

tle station to a flood-frequency curve or relation for that station, to 


3: Discussion open until September 1, 1959. To extend the closing date one month, 
2 written request must be filed with the Executive Secretary, ASCE. Paper 1994 is 
sart of the copyrighted Journal of the Hydraulics Division, Proceedings of the 
American Society of Civil Engineers, Vol. 85, No. HY 4, April, 1959. 
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determine the floods of various recurrence intervals. It was considered tha 

the use of mathematically fitted curves involving prior assumptions as to 

distribution might obscure the very relationships that were being sought, ani 

that graphically drawn curves, conforming as closely as possible to the 

original data at each station, would be most desirable. The plotting position 
n+l 


of each flood of record was computed by using the formula , where n is 


the number of years of record and m is the rank of the flood, starting with 1 
as the highest. This formula has theoretical basis and also agrees closely 
with the length of the period of record for the top flood, thus again conform- 
ing closely to the basic data. 

From the frequency curves, values were selected for the 1.2-, 2-, 5-, 10 
and 25-year floods at 170 stations, and for the 50-year flood at 40 stations. 
Where the length of record was less than 25 years, correlation with nearby 
stations was used to ascertain the appropriate plotting positions within a 25: 
year period for the annual peaks of record. The floods for each of the state 
recurrence intervals were then separately correlated with a large number ¢ 
physiographic and meteorologic factors. Initially, the relations were ana- 
lyzed by graphical multiple-correlation techniques combined with rank- 
correlation methods. The following parameters were studied: 


1. Drainage area 


2. Slope 
Channel slope No. 1 
Channel slope No. 2 
Average land slope 
Tributary channel slope 


wo 


. Water-supply index 
Mean annual runoff 
Mean annual precipitation 
Maximum 24-hour precipitation 
10-year 60-minute intensity 
Total rainfall volume in 6 large storms 


4. Shape factor 
L2/A, equivalent of basin length divided by basin width 
L/A, basin length divided by drainage area 
2al, the summation of small subdivisions of the drainage area, eac 
multiplied by the distance of travel to the gaging point. 


5. Mean elevation 
6. Percentage of lakes, swamps and reservoirs 


7. Stream density 
8. Temperature (January average) 


Most of the physical characteristics were computed as described by “3 
Langbein [1947]. “Channel slope No. 1” is one such characteristic. It is th 
mean slope of all parts of a channel draining at least 10 percent of the total 
drainage area,—thus it may include parts of the principal tributary stream: 
“Channel slope No. 2” is the same as that described by Bigwood and Thoma 
[1955]; the weighted main-channel slope is averaged with the tributary 
channel slopes as described by these authors. ; 
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he mean annual precipitation for each drainage area was Supplied by the 
». Weather Bureau as computed in the preparation of Hydrologic Atlas 

7 [Knox and Nordensen]. Other precipitation factors were based on vari- 
| Weather Bureau publications or on Geological Survey flood reports. 

y of the hydrologic factors, such as “channel slope No. 1” or “stream 

ity” were available from Water Supply Paper 968-C, Langbein, [1947], 
ugh not for all the 170 stations usable in New England. Because of the 

r involved in computing these, it was considered satisfactory, in the 
inary phases of the study, to compare the various factors and test 

r significance using only the 85 to 95 stations for which most of the 
rologic factors were already evaluated. 
in the graphical procedure used, Qn, the n-year flood, was first plotted 
log-log paper) against the drainage area. Through the plotted points, an 
rage curve was drawn, and the residuals from this curve were determined 
each station and ranked in order of magnitude. The values for all other 
ameters were then ranked and the differences between ranks were used 
ompute the Spearman rank correlation coefficients[ Wallis and Roberts, 
6]. The parameter with the highest correlation (considering only one 
nm each numbered group in the preceding table) was then used as the next 
iable in the graphical correlation. After the proper adjustments, a plot 
adjusted) first residuals versus the first parameter values led to a set of 
ondary residuals, which was in turn tested against the remaining parame- 
[Table 1 shows the results of graphical correlation for the 2-, 10-, and 
year floods for all significant parameters. This shows the decrease in 
standard error when successive independent variables were included. In 
h case, after 3 or 4 independent variables were used, further variables 
not appreciably decrease the standard error. It is notable that factors 
juating tributary slope alone were not found to be significant, once the 
ct of the main channel slope, as expressed in channel-slope factor no. 1, 
accounted for. The graphical standard errors were determined as the 
ye within which two-thirds of the plotted points lie. Actually these are 
standard errors of estimate, but are approximations adequate for these 
oratory studies. The standard errors are shown in terms of the original 
units and in the equivalent percentages these represent. One log unit 
ils one log cycle and a decimal part of a log unit is that same decimal 
of the linear distance representing one log cycle. 
.S would be expected, the standard error increases with the increase in 
size of flood.. It will be noted that for each size of flood, the standard 
)r with drainage area (A) drops sharply after channel slope (S) is in- 
ed, and is improved by only a small amount thereafter. The average 
entage improvement with S is 22, 38, and 19 percent respectively, 
reas the other parameters after A and S show combined improvements 
=xceeding 6 percent. 
- will be noted also that water supply indices, R, P, andI, show but small 
st. The relatively minor influence of rainfall is attributed to the uniform 
srn of precipitation throughout New England and to the unsatisfactory 
fall parameters so far available. 
he importance of channel slope is shown by the large percentage of the 
r which it accounts for. However there has been no unique or universal- 
‘cepted way of evaluating channel slope. Some hydrologists have used the 
drop from the head of the longest watercourse. Others have used 
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Table 1.--Results of Graphical Multiple Correlation 


Standard error 


Q Log pian = ila Bate 3 
E units e019 R- ee ee eee 
A 0.24 74 42 58 
eas .155 43 30 36 
A Sy) BR .15 Al 29 35 
| ae ee: fay sz 35 25 30 

210 
A ceul 104 51. 78 
Ors ay: 48 32 40 
Veale Pi Ee 42 29 36 

Qo5 
A .26 82 45 64 
A, Ss .19 55 35 45 
oe el 48 32 40 
hs LEG AR .165 47 aa 39 


A=Drainage area; S=Channel slope No. 1; R=Mean annual runoff: 
P=Mean annual precipitation; I=l0-year, 60-minute maximum 
precipitation; L/A=Shape factor. 


weighting methods which evaluate the slope all along the main channel. Sor 
have used parts of the main channel, such as the lower three-quarters. Sti 
others have in some manner combined the slopes of tributary streams with 
the main channel. Part of the difficulty is that methods of defining the mai 
channel or of differentiating the tributary streams are not fixed. 

Channel slope No. 1 was found to be the best slope factor of those evalu- 
ated. Its computation involves the determination of the points on the large: 
streams at which 10 percent of the total area is drained. This is a laborio 
task, particularly in a regional study where it may have to be done for 
hundreds of stations. It was therefore considered desirable to attempt to f 
a simpler, and possibly better channel slope factor. It was also considere 
desirable to separate the effects of tributary streams from those of the ma 
stream. It was decided (1) to use only the main channel in a channel-slope 
parameter, leaving the tributary slopes for separate consideration, (2) to ¢ 
fine the main channel, above each stream junction, as that channel draining 
the largest area, (3) to make an exhaustive study to find a way of expressir 
the main-channel slope that is most closely related to peak discharge. 

Topographic maps were used to develop channel profiles for all 170 sta 
tions in New England. At the upstream end each profile was extended to th 
drainage divide beyond the end of the stream shown on the topographic maj 
Distances were measured from the gaging point to each contour crossing ( 
cept where these were too dense) and channel profiles were plotted from 
these figures. 
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was considered that the most upstream part of the stream, in the steep 
jiwaters, might affect the slope out of proportion to the volame of water 
‘ished by the headwater area. At the downstream end the slope might not 
ndicative of that affecting the size of peak discharges because gaging 

8 are usually located at riffles or natural controls. It was therefore 
julated that the part of the main channel whose slope would best corre- 
/with peak discharge would be the one excluding the extreme headwater 
th, and possibly some of the extreme downstream end. 

The total lengths of the streams were divided into tenths below 0.7 of the 
ance from the gage and half-tenths beyond that. The object was to com- 
slopes for all possible portions of the channel, for example, 0.95 to 0.0, 
to 0.1, 0.6 to 0.4, etc., the figures representing in each case the fraction 
1e total distance upstream from the gage. There were 91 such combina- 
Ss. Two additional slope factors were computed. These were the con- 

its in the regression equation relating the logarithm of the rise to the 
-rithm of the distance from the gage, equivalent to drawing a straight line 
ugh the channel profile plotted to a logarithmic scale. These two factors 
resent an integrated slope for the entire main channel, rather than mere- 
1e Slope between two points. 

[o determine the optimum slope factor, multiple correlation analyses 

e to be made with the 1.2-, 2.33-, 5-, 10-, 25-, and 50-year floods as de- 
jent variables and drainage area and each of the 93 slope factors as the 
pendent variables. The change from the 2-year to the 2.33-year flood 
made at this time in order to use a value corresponding to the mean an- 
flood in common uSe; results are comparable. The best slope factor 

ld be defined at the point of minimum standard error. The correlation 
yses were such a prodigious task that the data were prepared for auto- 
ic solution by a digital computer, the Datatron 205. The programming 
sred computation of all the slopes, the logarithmic slope factors, the 
dard errors with drainage areas alone, and the standard errors and cor- 
tion coefficients with drainage and each of the 93 slope factors. Data for 
Stations was used in these correlation studies. 

tesults of the Datatron solutions are shown in Figures 1 and 2. These 
res show contours representing equal values of the standard error. On 
ire 1, for example, the standard error using the slope between points 0.5 
0.4 of the total distance above the gage is 0.193 log units. Figure 1 shows 
he standard errors on which the contours are based whereas Figure 2 

vs only the contours. Note the similar patterns shown by the controus 
2ach size of flood and the minimum in the upper left of each figure. A 
thting process reveals that the slope between points 85 and 10 percent 

'e the gage would satisfactorily give the minimum standard error for all 
is, with little accuracy lost for any particular size of flood. Standard 
rs using the two logarithmic slope factors were higher for all floods 
those using the 85 to 10 percent slope. 

‘able 2 shows the standard errors (in log units and average percent) for 

| size of flood, using drainage area and the “85-10” percent slope as in- 
ndent variables. The improvement in each can be seen when slope is in- 
uced. A regularity shows up in the percentage of the error explained by 
e, which is close to 35 percent in each case. The standard error for the 
ear flood seems inconsistently small and the contour pattern is some- 
erratic. The irregularities may rise because the 50-year flood is based 
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FIG. | VARIATION OF STANDARD ERROR WITH SLOPE FACTC 


on 40 items, whereas the others are based on 170 items. They may also be 
due to the less dependable data available for the 50-year floods. 

Comparison is shown with the graphical study (using channel slope No. 1 
at the same stages. Results by the datatron and graphical method are not 
strictly comparable, first because of the difference in methods, second be- 
cause 85 to 95 items were used in the graphical analysis and 170 items int 
datatron analysis. However, the final datatron results, using the “85-10” 
slope, indicate practically the same standard error as the channel slope 
No. 1 used in the graphical analysis. This is gratifying because the “85-1( 
slope is far simpler to compute. 

Figure 3 shows the typical regression equation with drainage areg and 
slope and a plot of the regression coefficients. The consistent variation in 
regression coefficients is evident. (The coefficient, b, for the drainage ar 
might well be 1.00 throughout). The multiple correlation coefficients for 
each of the regressions range between 0.93 and 0.97. 

Opportunity still exists for improvement in the standard errors shewn i 
Table 2, by the use of other parameters employed in the graphical study o1 
by the improvement of parameters such as precipitation indices. Further 
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FIGURE 2: VARIATION OF STANDARD ERROR WITH SLOPE FACTOR 
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Table 2.- Multiple Correlation with Drainage 
Area and Slope (85 to 10%) 


Standard error 


Datatron Graphical % Error 
explained 
Q log ave. ave. by slope 
ete units L % Bh ES) “4 
A CO. 222 53 = 34 
Aes .152 35 € 
Qo 33 
A eoek Be 58 on 
ras, ELSES 36 *36 
Qs 
A agen 61 - 38 
yee 164 38 = 
Q10 
A 264 65 78 40 
ANS 165 39 *40 
Qo5 
A 296 74 64 34 
A 8 207 49 *45 
850 
A «200 56 = 34 
Dees .158 37 bs 
Ave. 35 


*Graphical standard error computed using channel slope No. l. 


work is proceeding along these lines. This report is presented as a progress 
report with a definite finding on one phase of flood-frequency analysis, and 
is not intended to provide the final answer. 

The *85-10” percent main-channel slope found so effective for New 
England may not prove to be the best for other regions but it should serve 
as a guide in similar studies made elsewhere. 
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FIG. 3 VARIATION OF REGRESSION COEFFICIENTS WITH SIZE OF FLOOD 
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| STORM WATER DRAINAGE IN THE CHICAGO AREA 
| 


Horace P, Ramey,! M. ASCE 


SYNOPSIS 


This paper, read before the Hydraulic Division of the American Society of 
ivil Engineers on February 25, 1958, is a discussion of past and recent 
ooding conditions in the Chicago area and some recommendations for their 
nprovement. 

The opinions and recommendations stated herein are based on fifty years 
f investigation and experience with these problems. The recommendations 
iscussed are not meant to be specific, but rather to assist future students of 
ie subject in arriving at definite solutions. 


Floods have plagued mankind since the great Deluge of 2349 B.C. That was 
out 4000 years ago. One might think that the Chicago area, much of it ona 
ontinental Divide, with drainage both into Lake Michigan and into the 
ississippi Valley, would be comparatively free from floods, but such is not 
e case. 

To protect Lake Michigan, the source of its public water supply, from con- 
mination by sewage, Chicago decided many years ago to keep its drainage 
it of the lake. This eliminated the quickest outlet for storm drainage, and, 
cause of the topography, placed the main outlet for such drainage at the 
uthwesterly corner of the Chicago area. This, and the generally flat terrain 
the region, has produced many problems, some of which are yet to be 
lved. 

In the matter of flood runoff, the Chicago region is served by three 
reams: the Chicago River, which drains the northern area of the city itself 
far south as 87th Street; the Des Plaines River, which drains the westerly 
eas; and the Calumet River, which drains the entire Calumet region south 
87th Street. 


ite: Discussion open until September 1, 1959. To extend the closing date one month, 
a written request must be filed with the Executive Secretary, ASCE. Paper 1995 is 
part of the copyrighted Journal of the Hydraulics Division, Proceedings of the 
American Society of Civil Engineers, Vol. 85, No. HY 4, April, 1959. / 
Chf, Engr., The Metropolitan Sanitary Dist. of Greater Chicago, Chicago, 
Ml. 


11 


12 April, 1959 HY 4 


A knowledge of the topography of the Chicago region is essential to the 
study of its flooding; and a study of the past furnishes the best data for the 
solution of future problems. 


Topography, Drainage Areas 


The most important stream is, of course, the Chicago River, which extend 
due west from Lake Michigan, a little more than a mile, to the Forks. From 
this point the North Branch extends north by northwest, about twenty-nine 
miles, to its source in the Skokie Valley, west of Highland Park and Lake 
Forest, in Lake County, Dlinois. The lower six-mile reach of the North 
Branch of the Chicago River has been improved for navigation and flood flow, 
to Lawrence Avenue. Beyond that point the North Shore Channel of The Metre 
politan Sanitary District of Greater Chicago extends eight miles to a con- 
nection with Lake Michigan, at Wilmette. This channel (dry weather flow ca- 
pacity, 1000 cfs) passes the drainage of the North Shore towns, Evanston to 
Glencoe, into the North Branch. 

About three blocks north of Lawrence Avenue a dam was constructed to 
hold the North Branch of the Chicago River in its natural state, hydraulically. 
Above this dam, the North Branch extends westerly and then northerly about 
eight miles and then divides into three small streams, the: West Fork, the 
North Branch, and Skokie River. The Skokie Marsh covered a five-mile reac 
of this latter river, west of Winnetka and Glencoe, and, in its natural state, 
served to retain certain flood water. Much of this marsh has been drained 
and the area used for residential building; and lagoons have been constructed. 
to control to some extent the flow of water in the Skokie River. 

The South Branch of the Chicago River, extended from the Forks (Lake 
Street) due south about two miles, and then southwesterly two miles further, 
to Ashland Avenue and 26th Street, where it divided into the South Fork and 
the West Fork. The South Fork extended south one and one-half miles to 39tt 
Street and Racine Avenue. The West Fork extended west about two and one-— 
half miles to Kedzie Avenue and 32nd Street, the easterly limit of a slough — 
known as Mud Lake, which was the source of the South Branch. Mud Lake ex 
tended southwesterly, about six miles, to Harlem Avenue and 49th Street, 
where it connected with the Des Plaines River. 

The Main Drainage Canal was constructed from the West Fork of the Soutl 
Branch of the Chicago River, at Damen Avenue (one-half mile west of Ashlan 
Avenue) to extend southwesterly across the Continental Divide just south of — 
Mud Lake, past Summit, and down the Des Plaines Valley, to connect with the 
Des Plaines River at Lockport, Dlinois. ‘ 

co ae drainage areas of the Chicago River (above Summit) are 309 squa 
mules. 3 

The Des Plaines River, rising in southeastern Wisconsin, west of Racine, 
flows southward, parallel to and about ten miles west of the west shore of — 
Lake Michigan, to about 49th Street and Harlem Avenue, Chicago (one mile — 
south of Riverside), and then makes a sharp bend to the southwest and con- ~ 
tinues past Lockport and Joliet, to Dresden, where it unites with the Kankake 
River to form the Dlinois River. Drainage area of the Des Plaines River, 
above Riverside, is 643 square miles; above Lemont, 690 Square miles; and 
above Lockport, 713 square miles. _ 
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| The most important tributary of the Des Plaines River is Salt Creek, which 
pmes in from the west at about 39th Street. Salt Creek rises in Cook County, 
prtheast of Palatine, flows southward through the eastern section of DuPage 
jounty, past Bensenville, Villa Park, and Elmhurst, to Hinsdale; and then 
jows eastward past Western Springs and LaGrange Park, through Bookfield, 
nto the Des Plaines River. Its drainage area is about 152 square miles. 

) Flag Creek rises in the southern part of Hinsdale and Western Springs and 
Lows south to about 87th Street, then southwesterly and into the Des Plaines 
‘a about half-way between Willow Springs and Sag. Its drainage area is 

ut 20 square miles. 

The Little Calumet River rises in Laporte, Indiana, and flows westward 
nto Tllinois at about 170th Street, then northwesterly to a point in Blue Island 
bout one-half mile west of Ashland Avenue, at 136th’Street. Then it bends 
harply to the east and flows about seven miles to a junction with the Grand 
‘alumet River, at Torrence Avenue and 139th Street. From this junction the 
‘alumet River (drainage area 38 square miles) extends north about six miles, 
ast Lake Calumet, to Lake Michigan at 90th Street. From this same junction, 
he Grand Calumet River (drainage area 50.1 sq. mi.) extends eastward to the 
sake Michigan shore near the east line of Gary, Ind. Only a sand bar sepa- 
ates this river from the lake. The Grand Calumet is really a bayou, having 
nce been the outlet channel of the Little Calumet. 

The Indiana Harbor Ship Canal connects the Grand Calumet River with Lake 
fichigan, through East Chicago, Indiana. The drainage area of the Grand 
‘alumet, east of this canal, is 30.8 square miles; and west of it, 19.3 square 
tiles. 

The original watershed of the Little Calumet River, above Riverdale, was 
bout 670 square miles. This included about 75 square miles in the eastern 
alf of the Sag Valley, drained by Stony Creek, a tributary of the Little 
‘alumet, and seven square miles in the Blue Island area. It also included 
bout 370 square miles, in Indiana, in the upper reaches of the Little Calumet, 
hich, since 1925 have drained directly into Lake Michigan, near Dune Park, 
idiana, through Burns Ditch. Deducting these areas, total 452 square miles, 
rom the original 670 square miles, leaves the present watershed of the Little 
‘alumet, above Blue Island (its junction with the Calumet-Sag Channel) as 
pproximately 218 square miles. 

The drainage area of the Little Calumet River, between Blue Island and its 
inction with the Grand Calumet, is eight square miles, plus thirteen square 
liles tributary to the Calumet storm water pumping station; total thirty-one 
quare miles. 

The present drainage area of the Little Calumet River, as reversed, to 
ow into the Calumet-Sag Channel, at Blue Island, is 306 square miles. 

The drainage area of the Calumet-Sag Channel, east of the Continental 
ivide, is about eighty-two square miles; and west of this Divide, about twelve 
juare miles. 

The most important tributary of the Little Calumet River is Thorn Creek 
rainage area 115 square miles), which rises southwest of Chicago Heights 
id flows into the Little Calumet in the western part of Lansing, Il. One of 

e principal tributaries of Thorn Creek is Butterfield Creek (drainage area 
yout thirty square miles), which rises in the southwest corner of Cook 

ounty, west of Matteson, and flows through Olympia Fields, and then south of 
lossmoor and Homewood, to join Thorn Creek in Glenwood. An easterly 
ibutary of Thorn Creek is Deer Creek (drainage area about sixteen square 
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miles), which rises in Crete and flows northerly through East Chicago 
Heights, and then northwest along the Glenwood-Dyer Road, into Thorn Creek 
in Glenwood. 

The area south of 159th Street and west of the linois Central Railroad is 
drained by the Calumet Union Drainage Ditch (drainage area about twenty-fou1 
square miles), which runs easterly from Kedzie Avenue, along the line of 
163rd Street, and into the Little Calumet River in South Holland. 


Summary—General Topographic Characteristics 


The area discussed herein may generally be described as flat and rolling 
and although there are many ridges and divides in the area, they are not obvi- 
ously distinct. 

The existence of a large flood plain is evident in those areas containing 
ponds and sloughs as well as the areas bordering many of the shallow stream 
beds. 

Most of the streams in the area flow parallel to the lake shore as a result 
of the terminal morraines created by the receding glaciers. 


Changes in Stream Conditions 


The first change from the natural conditions of these streams, in regard 
to storm flow, occurred in 1848, when the Illinois and Michigan Canal was 
opened, from the South Branch of the Chicago River at Bridgeport (Ashland 
Avenue) through the Des Plaines Valley, to Lockport. About a year later, the 
Calumet Feeder of the I. & M. Canal went into service, from the Calumet 
River near Blue Island, through the Calumet-Sag Valley to the I. & M. Canal, — 
at Sag. 

Between 1867 and 1871, the Illinois and Michigan Canal was deepened, fron 
Bridgeport to Lockport, for a gravity flow of 1,000 cfs from the Chicago 
River, in lieu of about 250 cfs previously pumped. 

In 1871, the Ogden-Wentworth Ditch was dug, from the Des Plaines River, 
through Mud Lake, to the West Fork of the South Branch of the Chicago River 
The purpose was to drain Mud Lake, to enhance local real estate values. Thi 
ditch was 20 feet wide and its bottom was below the bed of the Des Plaines 
River. In the spring of 1872, the Des Plaines River floods greatly enlarged 
the new ditch, depositing the soil in the I. & M. Canal. 

Thereafter, more flood water than normal came from the Des Plaines 
River into the West Fork of the South Branch. Even the dry weather flow of. 
the Des Plaines River followed the same course and then flowed down the I. & 
M. Canal. This condition continued until June 1877, when the Ogden Dam was 
built. 

The Ogden Dam was a row of sheet piling, 50 feet long, with crest at ele- 
vation +11.8 CCD, which was about two feet lower than the natural banks of — 
the Des Plaines River, and a foot or two higher than the surface of Mud Lake. 
It was built in the arm of the Des Plaines River which led to Mud Lake, west. 
of Harlem Avenue, near 49th Street. This dam was rebuilt in 1885, of stone- 
filled timber cribs, with its crest at the same elevation. The new dam serve 
until 1894, when the Riverside Spillway, 397 feet long, at elevation +16.25 


CCD, was built on the east side of the Des Plaines, a short distance north of 
Ogden Dam. 
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The Riverside Spillway was closed in 1908 when the Willow Springs Spil 
way was built, to replace it, and discharge a like amount of flood water fro 
the Des Plaines River into the Drainage Canal, instead of into the South 
Branch of the Chicago River. The Willow Springs Spillway was closed in 
1955, and since that date Des Plaines River floods have not affected the flo 
in either the Chicago River or the Drainage Canal. 


Chicago Drainage Canal 


The greatest change in stream conditions in the Chicago region was the 
construction of the Chicago Drainage Canal, 1892-1900, from the West For 
of the South Branch of the Chicago River, at Damen Avenue, twenty-nine 
miles, to the Des Plaines River at Lockport, This was later extended, 19( 
1907, four miles, to connect with the Upper Basin of the Dlinois and Michi 
Canal and the Des Plaines River, in the northern part of Joliet. This chan 
was designed, 1892, with a normal hydraulic capacity to flow 10,000 cfs of 
water, to carry the runoff of the maximum storm expected in the Chicago : 
for the next thirty years. Its capacity, in the rock section, fifteen miles, f 
Willow Springs to Lockport, is 14,000 cfs. 

Construction of the Drainage Canal was accomplished by the improvem 
of the South Branch of the Chicago River, to capacity to flow 10,000 cfs of 
water at a velocity of less than 3 feet per second. 


North Shore Channel 


The North Shore Channel was constructed, 1908-1910, from Lake Michi 
at Wilmette, southwest about two miles to Emerson Street, Evanston, and 
then due south on the line of Kedzie Avenue, six miles, to Lawrence Avenu 
It joined the North Branch of the Chicago River about three blocks north o 
Lawrence Avenue. Its nominal flow capacity was 1000 cfs. 

Prior to the construction of the North Shore Channel, the marshy land 
around Emerson Street, and southward, drained south through a small bro 
into the Chicago River, south of Foster Avenue. The North Shore Channel 
added no significant drainage area to the North Branch of the Chicago Riv 
but it did speed the runoff of flood water. 


Calumet-Sag Channel 


The Calumet-Sag Channel was constructed, 1911-1922, from the Little 
Calumet River near Blue Island, sixteen miles through the Sag Valley, to: 
nect with the Main Drainage Canal, at Sag, about three miles upstream fre 
Lemont. The nominal flow capacity of this channel was 2000 cfs. Its dept 
was 20 feet and its minimum width 60 feet. This channel is now being 
widened, by the Federal Government, for navigation, to a bottom width of | 


feet, at a depth of 9 feet below water. This will increase its nominal hydr 
capacity to 5000 cfs. 


Summary Regarding Man Made Changes—General Effects 


The areas discussed in this paper are densely populated and highly ind 
alized. Some of the stream characteristics have been altered by dredging 
the addition of many man made structures. 
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n effect of these man made changes has been to change, in some cases, 
utlying borders of the drainage basins. This, together with the general 
raphic characteristics of gently sloping plains and flattish slopes, pre- 

§ a drainage problem which requires bold and imaginative solution. 
onsiderable progress has been made in the past with the creation of the 
ing open channel system, as well as the sewer program presently under 

It is, however, obvious that the tremendous population explosion and 

id rate of urbanization require immediate action to remedy a deteriorating 

pon. 

: 


oduction to Flood History 


onsiderable material describing early floods in the Chicago region is 

lable in the form of historic rhetoric. 

ome rainfall data as far back as 1875 is available for use. This data is, 

ver, not sufficiently complete to afford any scientific analysis of their ef- 
The following briefly summarizes the general impressions of those re- 

led floods and periods of intense rainfall. 


Early Floods 


‘he first record of a flood in the Chicago region was March 29, 1674, by 
sxplorer priest, Marquette. He and two other Frenchmen and a number 
idians, on their way from Green Bay to the principal town of the Dlinois 
ans, had spent the winter, after December 4, 1673, near Damen Avenue 
the West Fork of the South Branch of the Chicago River. On March 29, 

-, he was driven from his camp by high water coming through Mud Lake, 

1 the Des Plaines River, with the spring break-up of the ice. He proceed- 
he next day, up the South Branch, through Mud Lake, without any portage, 
then on down the Des Plaines River, which was in flood (Brown’s 

1inage Channel and Waterway,” pages 103-104). 


Flood of March 12, 1849 


most disastrous flood occurred in the Chicago River, March 12, 1849. 
recipitation records are available, but this flood followed the pattern of 
ious floods coming in from the Des Plaines River, through Mud Lake, into 
jouth Branch. An ice jam formed in the South Branch, south (upstream) 

e Madison Street bridge and raised the water and floating ice two to three 

The jam broke at 9 A.M., March 12, and the rapid current swept away 
y vessel moored below that point. The mass of ice and sailing vessels 
=k the Madison Street bridge, wrecked it, and added the bridge to the float- 
lebris. Next, the Randolph Street bridge was swept away, as were a 
er of canal boats and sailing vessels moored in the South Branch, north 
undolph Street. 
he mass of ruin then swept into the main river, where the banks were 
with steamers, schooners, brigs, canal boats, scows, etc. The Wells 
st bridge was open and the middle of the stream clear, but the floating 
; of vessels and wreckage extended from bank to bank and soon swept 
both ends of this bridge and mingled it with the general wreckage. The 
1r vessels crashed together, entangling booms, bowsprits, and rigging. 
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smaller vessels were crushed. The ice, in places, was wedged several 
above water level. 
(he Clark Street bridge was completely swept away and a great many of 
larger vessels became wedged together a short distance east of Clark 
pet, and with the general wreckage formed a dam, by 12 o’clock noon. 
jve this dam, ice and wreckage accumulated to considerable height; and 
pw it, numbers of vessels were swept into Lake Michigan, many of which 
re wrecked. 
he dam still held at 10:00 A.M. the following day, but most of the ice had 
sed under the vessels and the water level had fallen about two feet. The 
sis was over. Fears that the ice in the North Branch would break up and 
to the damage, were not realized. 
Me flood was general and caused considerable damage in the Des Plaines 
Illinois Valleys. 


Other Floods 


Other floods occurred in Chicago, due to the following extremely heavy 
ns: 


September 9, 1875 3.44 inches 
July 25-26, 1878 4.14 inches 
July 6-7, 1879 3.25 inches 
November 11-12, 1881 3.35 inches 
November 5-6, 1883 3.39 inches 
March 25-26, 1884 3.26 inches 
June 2-3, 1885 3.44 inches 


However, the records available are not sufficient for any study of these 
ds. 
mary of Runoff from Past Floods 


For convenience in study, the principal data regarding past floods in the 
cago area are listed in the tabulation following: 


Rainfall Time Runoff 
Inches Hours cis 
August 2-3, 1885 6.33 ad) ia 
May 6-8, 1892 ? ? 10,000 
August 11, 1908 4,34 24 9,000 
April 29, 1909 ese 36 8,900 
August 14-15, 1909 Shales 24 10,000 
March 17-18, 1919 (apes! 48 10, 600 
May 3-6, 1919 (remit Lhe 9,000 
Sept. 11-13, 1936 4,58 39 15,000 
July 1,5. 1938 4,14 16 12,000 


July 6, 1943 3.94, 54 24 16, 000 
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April 5-7, 1947 5.27 24 29,100 
March 19, 1948 3.50 13 24,500 
July 14-15, 1949 4.64 24 16,000 
June 2, 1950 4,38 24 13,000 
July 16-17, 1950 5.40 24 16,000 
May 10-11, 1951 3.91 24 16, 000 
October 9=11, 1954 6.66 28 27,000 
7.20 47 
July 12-13, 1957 6.00 6 28,000 
6.34 23 


The conclusion to be drawn from a study of these floods is that in the 
Chicago area within the past thirty years the runoff, from one cause or an- 
other, has increased at least two and one-half times. The Chicago Drainage 
Canal, planned in 1892 to carry the runoff from the greatest flood which mig 
occur in the Chicago area, then estimated at 10,000 cfs, has now become in- 
adequate. Many changes have taken place in the Chicago area since that tim 
sixty-five years ago, and the whole storm drainage problem must be re- 
evaluated. Apparently outlet capacity of 30,000 to 35,000 cfs will be needed. 


Floods in the Des Plaines River 


The Des Plaines River has a drainage area of 643 square miles, above 
Riverside. Of this, some 152 square miles are in the watershed of its prince 
pal tributary, Salt Creek. The drainage area of Salt Creek, above the Fuller 
burg dam, just north of Hinsdale, is about 122 square miles. | 

From the Hoffman dam, in Riverside, the Des Plaines River extends thirt 
one miles north to the Cook-Lake County line, just north of Wheeling, and or 
fifty-seven miles further to its source in Wisconsin, west of Kenosha and 
Racine. Downstream from Riverside this river extends southwesterly, abou 
twenty-six miles, to Lockport. The slope in the flood plain of the river, and 
in the river itself, averages about 1.0 feet per mile between Riverside and tl 
north line of Cook County, and about 1.9 feet per mile between Lockport and 
Riverside. 

The average cross section of the Des Plaines River, between Riverside a 
North Avenue, for dry weather flow, is slightly less than 400 square feet; an 
the section between North Avenue and Wheeling is probably less. This mear 
a channel about 150 feet wide and less than 3 feet deep. } 

The Forest Preserve District of Cook County, in conjunction with the Sta’ 
of Dlinois, has adopted a project to control the low water flow of the Des 
Plaines River, which at times is as low as 5 cfs, for weeks on end, in the 
thirty-one mile stretch between Riverside and the Cook-Lake County line, © 
This involved the construction of twelve to thirteen low head dams, to provic 
at least two feet of water above each dam. Dams at Dempster Street, Devon 
Avenue, and Touhy Avenue were built in 1946. The Hoffman dam, in Riversi 
was rebuilt in 1950-1956. Dams at Hintz Road, one mile south of Wheeling, 


and at Foundry Road, were to replace existing dams. One other dam was bu 
in 1957, ; 
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t least twenty-nine highway and nine railroad bridges cross this stretch 
e Des Plaines River, the average horizontal clearance of such bridges 

g about 150 feet. The average loss of head through these bridge openings 
imes of high water, may be 3 to 4 inches per bridge. The flood stage at 
th Avenue may therefore be 3 to 4 feet higher because of losses through 
e restricted openings. The flood stage at Wheeling may be increased 

€ as much. Accumulated drift at bridges and other narrow places has, no 
bt, caused further losses, for there is no record of this river having been 

ed during the past twenty-five years. 

Some recent floods in the Des Plaines River, as measured at Riverside, 
as follows: 


Maximum Runoff 
Riverside Discharge Cis 
Gauge cfs Devas deem. 

April 21, 1881 - 13, 500 21.0 
February 9, 1887 2022 10,324 16.1 
June 24, 1892 7a) os) 13, 000 20)..2 
March 25, 1904 21.45 10, 780 16.8 
March 7, 1908 21,47 10, 810 Wags: 
May 1, 1909 Plead 1,620 ise 
March 18, 1919 (ss eel 7,450 .? D6 
May 9, 1933 5,200 8.1 
pepts 17, 1936 5,160 8.1 
July 3, 1938 5, 350 $23 
March 16, 1944 4,600 et? 
April 6, 1947 3,780 aa) 
March 20, 1948 6,510 Mahe 
April 26, 1950 6, 340 Ciao 
October 11, 1954 6, 340 9.9 


[The measurement, April 21, 1881, was made by Lyman E. Cooley, of flow 
raweir. The measurement, February 9, 1887, was made by T. T. 

nson, using rod floats. : 
[he flood in April 1881 was due to the melting of ice and snow, amounting 
bout one foot of water, which had accumulated on ground which was practi- 
y saturated when winter set in. ; 

‘our floods in Salt Creek, as measured at the Fullersburg dam, in 1876, 
1, 1886, and 1887 produced flows of about 2800 cfs. In these cases, the 
ff was about 22.6 cfs per square mile, from the 122 square miles in the 
srshed above Fullersburg. The actual measurement in 1881 was 2760 cfs. 
the floods of 1881, 1892, and perhaps of 1909 may be classed as great 

ds. It is fortunate for the residents of the Des Plaines Valley that no 

it flood has occurred in the Des Plaines River since 1909. The highest 

sr levels in this time were in March 1919; and these levels were about 
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3.5 feet below the levels due to the 1881 flood. The recorded discharge at 
Riverside, March 18, 1919, was 7450 cfs. However, with a discharge of 651 
cfs March 20, 1948 and 6340 cfs April 26, 1950, the water levels in the Des 
Plaines River between Madison Street and North Avenue closely approached 
the high levels of March 1919. This indicates some obstruction to flow and 
the need for clearing the channel. 

A report on the “Improvements of the Des Plaines River” was made July 
1924, by Ernest L. Cooley, Hydraulic Engineer, to Edward J. Kelly, Chief 
Engineer of The Sanitary District of Chicago. For the critical reach betwee 
the Hoffman dam, at Riverside, and North Avenue, the improvement recom- 
mended was as follows: 


Lower the crest of the Hoffman dam from +25.56 to elevation +20.00 
C.C.Ds 

Excavate the channel from the Hoffman dam, 4800 feet to the junction 
with Salt Creek, to a bottom width of 300 feet and a bottom grade of 
+15.00; 

Excavate the channel, from the Salt Creek junction to a point about a 
mile north of North Avenue (Elmwood Park), to a bottom width of 200 
feet and a bottom grade rising to +17.00 at 22nd Street, +20.00 at Lake 
Street, and +22.00 at North Avenue; 

Construct a sub-dam at 22nd Street with crest at elevation +22.00; 

Lower the crest of the Public Service dam, below Lake Street, from 
+28.54 to elevation +25.00. 


This plan, estimated to cost several millions of dollars, would lower the 
flood water stage by 5 or 6 feet, from the Riverside dam to North Avenue ar 
beyond. The dams would cause low water boating pools, 18,000 feet long 
through Riverside to 22nd Street, 16,000 feet long from 22nd Street to Lake 
Street, and 15,000 feet long from Lake Street to Elmwood Park. 

This was a good plan in 1924 and it is a good plan today, with proper dl 
fications, including its extension to the north line of Cook County. Recent — 
studies indicate that the bottom grade of the channel should be about +13.00 
C.C.D. from the Hoffman dam to the mouth of Salt Creek and then rise, on a 
slope of 0.002, to about +17.20 at 12th Street, +21.20 at North Avenue, +33.0 
at Golf Road, +36.00 at Lake-Euclid Avenue, and +45.00 at the north line of 
Cook County. The depth, below high water, would be 10 feet. Bottom width 
the new channel would be 300 feet, from the Hoffman dam to the mouth of Sz 
Creek, then about 230 feet to North Avenue, then 210 feet to Higgins Road, 
diminishing to 150 feet at Palatine Road. Changes in channel widths would | 
made at points where tributary creeks enter the Des Plaines River. To ma 
tain desired low water levels, sub-dams should be constructed, at proper 
intervals. To lessen the rise in water level over these dams, to pass high 
flood flows, the dams could be built diagonally instead of directly across th 
stream, and thus be given greater crest lengths. Growth of shrubbery, whit 
would be obstructive to flood flow, should be restricted. 

It is believed that channel improvements, such as those suggested abov 
would enable the Des Plaines River to carry any flood which may Le 
be expected in the future. The depth of cut, below the bed of the present — 
channel, would probably average 4 to 5 feet, between the Hoffman dam and _ 
North avenue. and perhaps 3 feet between North Avenue and the County Lin 
The construction cost might be $6,000,000 to $9,000, 000. 
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provement would also be needed in Salt Creek, Addison Creek, Silver 
ek, and other tributaries. Reports have been made, by the Division of 
erways of the State of Illinois, indicating that such construction might 

t, as follows: 


| 


Addison and Silver Creeks $4, 790, 520 
| Flag Creek 1, 836, 170 
SHULS WEISS Sa Ot ao, eSNG 5,000, 000 


(Flag Creek (drainage area about 20 square miles), which rises in the 
ithern part of Hinsdale and Western Springs, north of 47th Street, is a 

all stream, about seven miles long. The Illinois Toll Road (now under con- 
ction) parallels Flag Creek, from 71st Street, north to about 47th Street, 

l crosses it. Considerable flooding has occurred, in the past, in Hinsdale 

1 Western Springs, because of the inadequate natural channel of Flag Creek. 
In recent years a considerable amount of building has been permitted on 
‘lands subject to overflow from high floods in the Des Plaines River. Inno- 
it purchasers of houses in such areas have suffered unexpected and un- 
rranted damage from high water. No public agency should approve such 
struction, regardless of the pressure applied by subdividers and builders. 


tle Calumet River 


A study of floods in the Little Calumet River indicates that storm runoff 
es in the Calumet area have at least doubled in the past fifty years. The 
inage area of this stream has been subject to several changes, over the 
irs. The original 670 square miles of the Little Calumet above Riverdale, 
luded seventy-five square miles in the Sag Valley and seven square miles 
und Blue Island, now tributary to the Calumet-Sag Channel, as well as 370 
are miles in Indiana which now drain to Lake Michigan through Burns 

ch. 

Records of the worst floods are as follows: 


Drainage 
Area Rate per 
Runoff Sq. Mi. Sq. Mi. 

Riverdale 

March 6-7, 1908 11,000 cfs 670 16.4 cfs 
North Harvey 
» June 11, 1939 3,440 cfs 198 17.4 cfs 

March 16, 1944 4,305 cfs 198 21.7 cfs 
Blue Island 

April 5-7, 1947 9,100 cfs 306 29.8 cfs 

March 19, 1948 5,000 cfs 306 16.3 cfs 

October 9-11, 1954 6,700 cfs 306 21.9 cfs 

July 12-13, 1957 7,200 cfs 306 23.5 cfs 


\ new lock and control gates are to be built in the Calumet River as part of 
project for widening the Calumet-Sag Channel, now under construction. 

$ will leave only the thirty-eight square miles of the Calumet River water- 
j tributary to Lake Michigan. The remaining 268 miles of the present 306 
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square miles of watershed of the Little Calumet, Grand Calumet, and Calum 
will be tributary to the Calumet-Sag Canal. 

For present flood relief into Lake Michigan, a channel 3, 500 to 4,000 
square feet in cross-section is available. In the future, only the capacity of 
the control gates at the new lock will be available. 

It would appear that the minimum capacity of the Calumet-Sag Channel, fe 
the passage of the runoff from a great flood, should be about 10,000 cfs. Thi 
much, or more, could come from a runoff of 30 cfs per square mile, from tk 
306 square mile watershed east of Blue Island and the 96 square miles in th 
Calumet-Sag valley. About 20 per cent of the 221 square miles within the 
city limits of Chicago, or 44 square miles, south of 87th Street, are includet 
in the area draining to the Little Calumet River and the Calumet-Sag Channe 

The channel of the Little Calumet River, upstream from Blue Island, is 
quite inadequate for storm water flow. In 1941 the State of Illinois improvec 
this channel for a 50-foot width, from the Indiana State line to the mouth of 
Thorn Creek; and in 1946 extended this 50-foot improvement to Roll Street, 
in Blue Island. Under present conditions of greatly increased runoff due to 
the considerable building in the area, it appears that the width of this channe 
should be more than doubled and perhaps the depth should be increased. Op 
ings at the many bridges should be greatly increased. A pending improvemée 
in Thorn Creek, already authorized by the State of Dlinois, will increase the 
rate of storm water flow; and direct attention to the need for further im- 
provement in the Little Calumet River above Blue Island. 


Summary—Floods: Need for Remedial Measures Obvious 


The foregoing has outlined chronologically, the history of storms and re- 
sultant floods in the Chicago region. As can be seen, there exists some dat: 
of recent years which can be used to appraise and investigate a few past 
storms and their relation to flooding. It becomes more obvious upon such é€ 
amination that immediate remedies must be found for the flooding situation 
the area. The flood frequency is becoming greater with successively small 
rainfall and promises to become still greater with further urbanization and 
growth, 


City of Chicago 


The City of Chicago embarked on a program of construction of storm rel 
sewers, in 1947, expected to entail an ultimate expenditure of $165,000 ,000, 
to catch up with existing needs. Three bond issues to provide funds for this 
construction have already been approved by the voters, as follows: 


L942Cicatee te lek Bieta = Boe $58, 160, 000 
LOS Sees. Shee et ee ee 30,000, 000 
Ee oS! eee rile Sk ee ee eee ee 8,000, 000 
Lotalvapprovedsrens se .opsne - - - $96,160,000 
PROtUPe-Geeos cliente tna oe oe 0 OOOnO00 


Total Pip «iahaundl co\e bs apse SIL6 S1GOS000 
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is hoped that this work can be carried on at the rate of about $15,000,000 
ear, 1957 and after. 
n addition to the work by the City of Chicago, The Sanitary District of 
ago expended $47,367,000 on its system of South Side relief sewers (to 
pmber 31, 1956); and is currently spending about $5,200,000 on an ex- 
tion of the Racine Avenue Pumping Station, for pumping storm water only. 
|total expenditure, by the Sanitary District, for this storm drainage will 
ibout $52,570,000. 
his system of sewers with the existing 39th Street Conduit was designed 
rry 5500 cfs of storm runoff, from 27.6 square miles, between 31st 
et and 87th Street, and east of Racine Avenue to Lake Michigan; or at a 
off rate of 200 cfs per square mile. These sewers have capacity to carry 
cis of sewage, in addition to the 5500 cfs of storm flow, total 5900 cfs. 
“he sewer outlet capacity of the City of Chicago, past and prospective, in- 
ling the South Side sewer system of the Sanitary District, is indicated in 
tabulation following. Runoff rates are based on the present 221 square 
2s of area within the City of Chicago. 


Runoff 

Years Sewer Outlet Capacity cfs per sq. mi. 
1900-1910 10,000 cfs ASS 

1920 14,000 cfs 64.0 

1930 20,000 cfs 91.0 

1940 30,000 cfs 13720. 

1950 35,000 cfs 159.0 

1957 37,000 cfs Pas0 

1960 41,000 cfs 187.0 

1970 55,000 cfs ; 250.0 


ewer outlet capacity, thirty years ago, was about ten times the dry 

her flow. The outlet sewers now being constructed are designed with ca- 

ties thirty to fifty times the dry weather flow. This will be reduced by 
water levels in the Chicago River and the Drainage Canals; but the 

age capacities in these large relief sewers will be considerable. 

hese present and prospective sewers are designed to flow their full ca- 

ties into the rivers or canals at a water level of +1.50 CCD. This water 

| may reach +4.00 in the Chicago and Calumet Rivers and the effect of 

would be to reduce the flow capacities of these sewers to 70 or 75 per 
of their designed capacities; or to reduce the total (1970) capacity from 

00 cfs to about 40,000 cfs. 

is not expected that the flow capacities of all these outlet sewers could 
be developed at the same time. The times of concentrations of floods in 
arious sewers and sewer systems range from less than thirty minutes, 
)three hours. However, a 6-inch rain in six hours, such as fell in the 

e Chicago area on July 12, 1957, could, through this completed sewer 
2m, deliver an amazing amount of storm runoff into the Chicago and 
met Rivers and the connecting canals. In July 1957 a lot of this went into 
asement of householders, because of the inadequacy of the local sewers. 
9 doubt these sewers, in particular, the recent auxiliary outlet sewers, 
been designed on a liberal basis. However, if they have not been greatly 
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over-designed, a study of how much flood water they can deliver to the outl 
streams may be a more realistic method of estimating the volume of a futu 
great flood, than the method outlined on previous pages, of estimating it fre 
the amounts of water which have been discharged through the Drainage Can 
in times of past floods. 

A study of these outlet sewers, in only the 221 square miles within the ¢ 
limits of Chicago, indicates the following: 


Percentage of Sewer Capacity - cfi 

Sewer Capacity Total 75% 5 

North of Lawrence Avenue 10% 5,500 4,120 2a 

Lawrence Avenue to Lake Street 17% 9,350 7,010 4, 

Total North Branch 27% 14, 850 TrPIS0 % 

South Branch Chicago River 30% 16,500 12,370 8) 

Total Chicago River 57% 31,350 23,500 Tay 
Drainage Canal, Western Ave. to 

Harlem Ave. 16% 8,800 6, 600 4 

Total to Main Channel 73% 40, 150 30,100 . 2am 

Calumet Area, South of 87th Street 27% 14, 850 115930 _% 

Total, City of Chicago 100% 55,000 41,230 27, 


Storm Flow—Computed from Sewer Capacity 


It may be assumed, as quite probable, that at least 50 per cent of the to 
flow capacity of these sewers will be developed simultaneously, in some fu 
flood, and that the storm runoff from the 221 square miles within the presi 
Chicago city limits will be at least 27,500 cfs. To this must be added the 
storm runoff from about 100 square miles, north of Chicago, in the upper 
reaches of the North Branch and of the North Shore Channe, including | 
Evanston and the other North Shore towns; and the runoff from about 260 — 
square miles tributary to the Little Calumet River, south of the city. Sucl 
runoff may reasonably be calculated at 30 cfs per square mile; and would 
amount to about 3000 cfs on the North Side, and 7800 cfs in the Calumet 
region. 
The 3000 cfs on the North Side, added to the 20,070 cfs total flow to the 
Main Channel, from the table above, would indicate a storm flow of at leas 
23,070 cfs in the Main Drainage Canal, above Summit. This runoff, from t 
308.9 square miles of drainage area above Summit, would be at an average 
rate of about 75 cfs per square mile. ; 

The 7800 cfs in the Calumet region, added to the 7430 cfs flow to the 
Calumet-Sage Channel, from the table above, would indicate a storm flow 
about 15,230 cfs in the Calumet-Sag Channel. This runoff, from the 306 si 
miles tributary to the Calumet-Sag Channel, would be at an average rate c 
about 50 cfs per square mile. Storm flow from the ninety-six square mile 
in the Sag Valley is ignored in this calculation, since much of this could hi: 
already passed downstream before the peak concentration in the entire : 
Calumet region. ; 

The total storm flow, under this method of calculation, would be 23,070 
from the area north of 87th Street plus 15,230 cfs from the Calumet regio: 
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8,300 cis. This can be compared with the conclusion of 30,000 to 35,000 
Eerived (page 25) from a study of the runoff from past floods. 


} 


roduction to Recommendations 


'The following discussions and tabulation comprise a program of action 

ired to increase the outlet capacities for the expected flood waters in the 

irs to come. They are by no means the total solution to the problem but 
make available sufficient outlet capacities for the rapidly expanding 

ver programs. It should be understood that these recommendations are 

sed on the data and experience taken from the operation of the open channel 

stem and are not meant to solve localized problems due to inadequate sewer 

. This program is suggested as that part of the total program which 

I fulfill the outlet capacity needs for flood waters for the next fifty years. 


Main Outlet Channels 


Since storm runoff in the Chicago and Calumet areas has approached 30,000 
on three different occasions in the past ten years (March, 1948, October 
44, July 1957), and since building and paving are increasing at very sub- 
ntial rates, it is apparent that the runoff from a great storm may be 30,000 
35,000 cfs, in the not distant future. It is obvious that outlet capacity for 
sh storm flow should be at least 30,000 cfs. 

Chicago is the one city, on the Great Lakes, in the fortunate position of not 
jing to drink its own sewage, or pass it on for others to drink. However 
ute it might be, the thought is repugnant. For the past thirty-five years, or 
ce August 1922, when the Calumet-Sag was opened, Chicago has kept all of 
sewage out of Lake Michigan, the source of municipal water supply. It has 
been possible to keep the peaks of flood flows from the Calumet Rivers 
of the lake, whenever such flood flows have exceeded the hydraulic capaci- 
of the Calumet-Sag Channel. The same is true of the North Shore Channel, 
ere storm runoff in recent years has had a great increase, and where the 
isequent water level has, on many occasions, been high enough to cause 

w over the restraining walls at Wilmette. On only two occasions (October 
44 and July 1957) has it been necessary to permit the Chicago River to flow 
) Lake Michigan, to prevent damage to Loop property, from high water. 

2se were on occasions of what might be called great storms. 

It may be that this policy of Chicago, cherished over the years, of keeping 
drainage out of its drinking water, may have to be abandoned; at least dur- 
great storms. It is recommended, however, that sufficient outlet capacity 
y be provided for storm flows so that relief discharge into Lake Michigan 

| be limited to the peak runoff from great storms. 

To that end, a hydraulic study has been made of the present outlet channels 
of how much their flow capacities should be increased, to meet the needs 
he near future, of discharge capacity of 30,000 cfs, assuming that any flood 
y in excess of that amount could be spilled into Lake Michigan. This study, 
ch should be considered as preliminary only, indicates a need for the 
owing improvements: 


Main Outlet, Drainage Canal, from Sag to Lockport— 

_ Widen present channel, from 160 feet to 300 feet, or more. 
This channel would carry about 30,000 cfs, with water surface -3.0 
at Sag and -8.0 at Lockport. 
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Calumet Outlet, Calumet-Sag Channel, from Blue Island to Sag— 
Present widening to 225 feet, by Government. This channel should 
carry about 10,000 cfs, with water surface -3.0 at Sag and +4.0 in 
Little Calumet River. 


Chicago Outlet 
Drainage Canal, from Willow Springs to Sag— 


Widen present channel, from 160 feet to 230 feet. This channel 
should carry about 20,000 cfs, with water surface -1.5 at Willow 
Springs and -3.0 at Sag. 


Drainage Canal, from Summit to Willow Springs— 
No improvement indicated. Present channel should carry about 
20,000 cfs, with water surface +0.1 at Summit and -1.5 at Willow 
Springs. 


Drainage Canal, From Western Avenue to Summit— 
Widen present channel 70 feet. The channel should then carry 
about 20,000 cfs, with water surface +1.5 at Western Avenue and 
+0.1 at Summit. Perhaps this reach should be improved to the 
standard dock section, with 240 ft. width between docks; and 
proper deepenings. 


Chicago River and Drainage Canal from Lake Street to Western Avenue- 
No improvement indicated. Present channel should carry about 
20,000 cfs, with water surface +4.0 at Lake Street and +1.5 at 
Western Avenue. 


These estimates are based on a maximum stage of +4.0 CCD in the Chic 
and Little Calumet Rivers, of -3.0 at Sag, and a minimum of -8.00 at Lock- 
port. Holding the water level at Sag at about -3.00 is very important. In th 
March 1948 flood, this level was about +1.60. In the October 1954 flood, it 
was probably much higher, but the gauge failed. In the July 1957 flood, it w 
+0.40. These high water levels at Sag cause an unreasonable decrease in tt 
flow capacity of the Drainage Canal above Sag; and in the Chicago River. — 

Further improvement in the Chicago River and in the Drainage Canal be: 
tween Damen Avenue and Summit could be easily accomplished by dredyia 
the central 100 feet of the cross-section to a greater depth. 

The improvement of these channels, as indicated, would involve the exe 
vation of approximately 7,200,000 cubic yards of Pock and about 4,300,000 
cubic yards of glacial drift, or earth, at a cost of about $30,000, 000. * mo 
detailed study would result in a better hydraulic balance between the va 
channels and sections of channels. 


Efforts 1955 to 1957 
: 
Since 1955, when the storm water drainage program within The Sanitary 
District of Chicago, outside the City of Chicago, then estimated to cost 
$72,000,000, was submitted to the Dlinois General Assembly, little has 
happened to panee the picture. 
The Metropolitan Sanitary District of Greater Chicago has caused the 
dredging of 117,325 cubic yards of material from the North Branch of hey 
Chicago River and the North Shore Channel, at a cost of $378,960. This in 
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nlarged this channel, but such maintenance dredging did restore the 
nal cross-section and improved conditions for flood flow. 
urrently, maintenance dredging of the North Shore Channel, from Foster 
jue to Wilmette, is under way. Some 250,000 cubic yards of silt and fill 
10 be removed, during the year 1958, at a cost of about $1,110,000, to 
ore the balance of the North Shore Channel to its original cross-section. 
bris and silt and some shrubbery were removed from the channel of the 
(Plaines River between the old Lockport Controlling Works and Ninth 
et bridge, in Lockport, and the original channel restored in 1956. Similar 
- is now under way between the Ninth Street bridge and the Elgin, Joliet, 
stor Railway bridge. This will be finished in 1958. 
ree additional sluice gates at the old Lockport Controlling Works were 
oped with operating machinery, 1956-1957, thus adding about 7,500 cis to 
lischarge capacity from the Main Drainage Canal, at Lockport. Current- 
tudies are being made of methods to obtain 10,000 to 12,000 cfs more out- 
‘capacity, probably at the Lockport Power House. Outflow capacity was 
2what inadequate during the October 1954 flood, when two of the sluice 
5 were blocked by some loose barges. 
he channel of Deep Run is being enlarged between Sixteenth Street bridge, 
»ecKport, and its junction with the Main Drainage Canal below the Lockport 
=r House. The improved channel will have a storm flow capacity of 6500 
This work will be finished in 1958, at a cost of about $500,000, and 
ld alleviate the worst flood conditions in Lockport, Dlinois. 
he City of Chicago continued construction on its long range program of 
liary Outlet Sewers, presumably at a rate of about $15,000,000 per year. 
nd issue of $8,000,000, for this work, was approved by the voters, in 
, leaving about $70,000,000 of the future financing yet to be approved. 


Estimates of Cost of Future Work 


) little has been done toward the alleviation of flood conditions in Cook 

ty, outside the City of Chicago, that the estimates of costs, submitted in 

, may well be repeated, with minor modification. One important and cost- 
1m, however, must be added, namely, the enlargement of the Main Outlet 
nel, at $30,000,000. 

he entire future storm drainage program, not yet financed, is outlined in 
vbulation following. 


rm Water Drainage Program - February 1958: 


Enlargement of Main Drainage Canal. ...+++++-+- $30,000, 000 
Improvement of Des Plaines River, from 

Riverside to north line of Cook County ....+++.+-> 9,000, 000 
Improvement of North Branch of Chicago River, 

north of Lawrence Avenue ...+++++*eee e228 4,000, 000 
Enlargement of North Ghore Ghannel. ai. 2 2-6 @ « e's 2,000, 000 


Enlargement of Stony Creek, and construction 


of tributary ditches and sewers... +++ +++ sees 8,000, 000 
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Improvement Calumet Union Drainage Ditch. ....+-.- 2,000, 000 
Improvement Feehanville and Weller Ditches......-.-- 1,000, 000 
Improvement Thorn Creek... . 1+ 2+ + +++ ee ee ee 1,000, 000 
Improvements along Addison and Silver Creeks ......- 4,000, 000 


Storm water outlets for Posen, Midlothian, 
Markham, ‘South: Blue Islands) Ja a ss ss ee 4,000, 000 


Storm water outlets, etc. for villages in 


Grand Calumet and Little Calumet areas ......2.2--s 6,000, 000 
Total Cook County, outside Chicago....... $72,000, 000 
Total City of Chicago, Sewers . 1... 6 «+s ae 70,000, 000 

Total storm water outlets....... $142,000, 000 


Storm water sewer systems for the various villages in Cook County, esti 
mated 1955 to cost $31,000,000, are not included in the above estimate, Th 
improvements of Salt Creek, in DuPage County, estimated 1955 at $5,000,01 
is of some importance to the Chicago area, but is omitted. 


CONCLUSIONS 


From a study of floods in the Chicago area, it is obvious that the presen 
outlet channels are not adequate to handle the storm runoff during a great 
flood. No change has been made in the Main Outlet, that is, in the Drainage 
Canal between Sag and Lockport, since 1900. In that time, the storm runoff 
from the Chicago and Calumet areas has probably more than tripled. Loca 
drainage has been and is being improved in many parts of the area, but thes 
improvements can not be fully effective until the Main Outlet is made adequ 

Therefore, the enlargement of the Drainage Canal might be considered a 
of first importance. This will be particularly true of the stretch of this cal 
downstream from Sag, since the widening of the Calumet-Sag Channel will, 
the next four years, cause a tremendous increase in storm flow from the © 
Calumet region. Enlargement of the Drainage Canal is estimated to cost — 
$30,000,000, or about the amount of damage estimated from one great floo¢ 
in October 1954. Damage from the flood of July 1957 was less, but a figure 
$10,000,000 has been mentioned. It is probable that the cost of flood pro- 
tection would be less than the cost of flood damage. 

All the projects listed in the program above, except the enlargement of 1 
Drainage Canal, are urgently needed for protection against the ordinary flo 
or that due to a ten-year storm. Residents of some of the communities in- 
volved were in pitiful distress during the great floods of 1954 and 1957. 

It is not the purpose of this paper to suggest specific improvements for 
handling storm water flow in any part of the Chicago area; but rather to ou 
line the problem and furnish background material and factual data which wi 
assist future students of the subject in arriving at definite recommendation 
Proper solution of the problem will require a construction program far bey 
the imagination of everybody except a few who have been close enough to it 
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3 its magnitude. Considerable time, at least ten years, will be required 
e physical construction, in an orderly fashion, of the needed works, 

> any method of financing which will not be too burdensome, It is a 

cam which must necessarily be carried out on an area-wide basis. 

nicago is a great city, although it is a young city, only 120 years old. 

aps it has grown too fast and some of its problems have been neglected. 

a become the greatest city in the world; and the solution of the storm 

tage problem, now, would enhance that prospect. If the powers that be 

race that essential issue, and act, future generations will bless them. If 
rocrastinate and postpone such really necessary improvements, future 

rations will condemn them for lack of courage. It is well known that inde- 

m is the greatest failing of the average office holder; but there is no time 

ndecision in a matter so important to a large and populous community as 

revention of flooding. The longer such matters are deferred, in a grow- 

ommunity, the greater the cost will be. ; 
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HYDRAULIC ANALYSIS OF SURGE TANKS BY DIGITAL COMPUTER 
: Nicholas L. Barbarossa,! M. ASCE 
: 


SYNOPSIS 


It is shown how certain surge problems have been successfully adapted for 
ution by an IBM 650 computer. Machine results for the Oahe surge system 
presented and compared with those obtained by independent manual 

thods of arithmetic integration. Machine-computation time for each of the 
es discussed averaged between one to one and one-half hours. Sufficient 

a are given to permit individual checks on machine analysis and compu- 
ons in a relatively short time. 

Computer results are discussed to indicate that treatment of the total 
blem of transients as one of separable parts in the analysis of surge-tank 
oility by conventional methods leads to inconclusive results. It appears 

- a comprehensive study by digital computer of the total problem of 

isients is feasible, but a machine programming job of great scope will be 
uired before fruitful results are available. 


INTRODUCTION 


ligh-speed, automatic digital computers have the potentiality of opening up 
horizons of understanding of problem phenomena. Machine solutions will 
n provide a broader perspective of problem phenomena than can be feasi- 
obtained by manual computations because of the relative rapidity with 

eh the effects of various factors and assumptions involved in a problem 

be tested and explored. By demonstrating how this advantage of a digital 
puter has been actually utilized for solving surge problems in a particular 
2, additional evidence is provided of the growing usefulness of computers 
ngineering work. 


: Di date one month 
; Discussion open until September 1, 1959. To extend the closing A 
written request must be filed with the Executive Secretary, ASCE. Paper 1996 is 
art of the copyrighted Journal of the Hydraulics Division, Proceedings of the 
merican Society of Civil Engineers, Vol. 85, No. HY 4, April, 1959. 

hf., Hydrs. Section, U. S. Army Engr. Div., Missouri River, Omaha, 
lebr. es 
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The study discussed herein is concerned with mass-oscillation phenomena 
in a surge system, which includes a pressure conduit and two restricted- 
orifice type surge tanks in tandem independently connected to the main con- 
duit (Fig. 1). The system supplies water to a reaction-type turbine. Water- 
hammer effects are not included in the study. Thus, the study is primarily 
concerned with the solution of surge-tank problems, involving load rejection, 
load demand, and hydraulic stability of the tanks under small load changes at 
minimum head. 

Given the initial-flow conditions and all other pertinent data, the values of 
various discharges and heads as a function of time for the unsteady flow con= 
ditions within the surge system, resulting from power-load changes, may be 
determined by conventional methods of arithmetic integration.(1,2)* Appli- 
cation of these methods to surge problems is based upon the fundamental as- 
sumption that the compressibility of water and the elasticity of flow bounda- 
ries may be ignored. This means that pressure changes in the surge system 
are transmitted with infinitely high velocities and that water columns in diffe 
ent parts of the system may be regarded as solids in evaluating positive or 
negative accelerative effects. The same fundamental assumption is used as ¢ 
basis for development of machine methods. However, since a digital computé 
is to be employed to obtain solutions, there is no need to make as many other 
simplifying assumptions as are ordinarily dictated by the time limitations of 
manual solutions. For example, acceleration heads within the risers, vari- 
ation of friction factors with changing Reynolds Number, tee losses, and fine! 
time increments in the step-by-step process may all be introduced with only 
comparatively small increase in the machine time required for problem so- 
lutions. . 

A list of the physical relationships required for machine solution of surge 
problems associated with the system depicted on Fig. 1 and the sign conventi 
that has been used are given in Appendix III. 


7 
Basis of Analysis and Computation 


Head and Discharge Relationships 


Based upon the fundamental assumption of instantaneous transmission of 
pressure changes, application of the energy principle to the upstream conduit 
(Fig. 1) for any instant of time yields** 


i 

hora Beet uly | 

ay R rE Ly i 

The velocity head for the upstream conduit is included in hy,, for conveniene 
also included are all form and surface-resistance losses between the a 


i 


voir and point (1). Relationships similar to Eq. (1) for other parts of the su 
system are given in Appendix III. 


All acceleration heads are defined by the conventional, one-dimensional 
formula 


heh 
aenek oat 


*The numbers in parentheses indicate references listed in Appendix II. 
**The symbols used are defined in Appendix I. 
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The sign of h, will depend upon the sign of the derivative dt 
By appropriate use of Eq. (2), Eq. (1) may be written in the form 


ag 
wiwe BA (p=, - by) (3) 


dt L uf L, 
Similar relationships may be developed for other Q’s. The diameter of the © 
main conduit has the same value for the lengths Lj, L3, and L5. All corre- 
sponding dimensions of the upstream and downstream risers and surge tanks 
are the same. 
The net head on the turbine is obtained from the relationship, 
2 
Vv 
= ee (4 
H ite + Seale Ty 
In Eq. (4), the comparatively small velocity-head lost in the tail-race is ig- 
nored. Moreover, Tw is assumed to be constant during any one solution. 
Application of the continuity principle at each tee requires that 


Q) FO = (5 


at each instant, even when flow reversals take place. Differentiation of Eq. 
(5) gives 


ETRE Re fl ers (7 


A similar expression may be obtained by differentiating Eq. (6). Equations o 
the form indicated by (7) are useful for eliminating the derivatives from 
equations of the form indicated by (3). Actual use is shown in Appendix III. 

The relationships between the water surfaces in the surge tanks and the 
instantaneous flows in or out of the tanks are obtained from the consideratior 
that inflow minus outflow equals the change in storage during dt: 


a(Ws, ) Q, 


= (8 

dt Ap 
a(Ws.5) Q, (9 
Qil ete ake 


For cases when the water surface in either tank drops below the level Eg, Aj 
in either Eq. (8) or (9) is replaced by Ay. 


Head Losses 


All surface-resistance losses hg in feet are evaluated by the Darcy- 
Weisbach equation, 


L 
ne = 2 5 LG | (10 

2gA 
This equation is applied not only to the main conduit but also to the risers. 
Surface-friction losses within the surge tanks are ignored because, ordinaril 
they are negligibly small. 
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\In some installations a part of the total length L, of the upstream conduit 

1 be lined with concrete and the remainder with steel. Provision for this 
been made in the machine program. Let Le = length that is lined with 

ine and Lg = length that is lined with steel, so that Lj = Le + Lg. Desig- 


ing fe and fg as the friction factors for the concrete and steel lengths, re- 
etively, and utilizing Eq. (10), the total surface-resistance losses from the 
ervoir to point (1) are given by 


| hp = (f.L4 21 ) al (11) 

1 S$ 2gDA 
Except for the tee losses, all form (or shape) losses hg in feet are evaluat- 
by an equation of the type, 


lala 
h — K eee 
s De® (12) 


which K is assumed to be constant for a particular flow-boundary configu- 
tion. Eq. (12) indicates that hg will have the same sign as Q. The total 
ad loss* due to various boundary forms from the reservoir to point (1) is 


ren by 


Q, 12 
on = (K+ Kot K+ K+ Kt Kt eS ee) wild | — (13) 
2g 
‘ictly, all K’s in Eq. (13) will have one set of values when Q is positive and 
il have a different set of values when Q, is negative (flow toward the reser- 
ir). In this study the K’s are assigned one set of values and are assumed 
remain constant during any one machine run. (If deemed necessary, the 
chine program can be modified to allow the K’s to vary with the direction 
Q, and also with the magnitude of the velocity head.) 
Combining Eqs. (11) and (13), the total head loss from reservoir to point 
is 


He =[q+ a! (fh, £2.) [24] % (14) 


Head losses (other than tee losses) for other parts of the main conduit are 
ailarly determined, using analogous notation (Appendix III). 
In evaluation of the losses in the risers and at points (5) and (12), attention 
st be given to the direction of Qg and Q4 and to whether or not the water 
‘faces Ws, and Wsg are above or below the elevation E, at any instant under 
sideration. The total head loss for the upstream surge tank and riser 
stem is 

b, = (C+ Ktot a ) [a] 4 (15) 

5 . 2 eee 5 

The first term in Eq. (15) evaluates: (a) the loss of velocity head when the 
er is delivering water to the surge tank; (b) the velocity head in the riser 
| the entrance loss when the surge tank is delivering water to the riser. 
> second term in Eq. (15) evaluates the surface-friction loss in the riser. 


rictly, hg is more than a head loss, since it includes the velocity head in 
e upstream conduit. K, is equal to one in all cases. 


“ 
k 
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Eqs. (I-27) and (I-28) in Appendix III indicate how Cy and Kp are evaluate 
In the two expressions for Kg (Eq. (III-28)) allowance is made as to whether 
or not the total friction boundary has a length R or something less than R. 
The third term in Eq. (15) determines the head drop across the orifice and : 
based upon the standard-orifice formula. In this study Cp is taken as being 
constant during any one machine run. Actually, even if the orifice may be 
considered a thin-plate orifice as such, it is likely that Co will vary as Q2 
changes owing to the condition that the orifice plate is probably not sufficier 
distant from major flow disturbances in most installations. However, Co m 
be easily introduced into the machine program as a function of appropriate 
variables, if desired and if data are available for realistic evaluation. Othe 
wise, Cy may be varied from one run to another to see what effect different 
values have on the hydraulic transients, as has been done in the Oahe study. 
Head losses hy,, in the downstream riser and surge tank are obtained in the 
same manner as used to obtain hi,9- 


Tee Losses 


In determining the tee losses at any instant of time, the three Q’s at eacl 
tee must be known in both magnitude and direction at that instant. Evaluatir 
of the tee loss will then depend upon which of the possible flow patterns indi 
cated on Fig. 2 applies at that instant. Provisions have been made in the 
machine program for evaluating the tee losses at each instant of time for a 
flow patterns shown on Fig. 2 except the last two (i and j). (When flow 
patterns 2(i) and 2(j) occur, the tee loss is assumed to be zero.) For these 
two cases, it is uncertain that appropriate data for evaluation exist. Should 
data be available, the program could be readily modified to accommodate 
these two cases. Data are available(4) to evaluate cases (a) through (h) on 
Fig. 2. Obviously, cases (e) through (h) are special cases of (a) through (d) 

The tee losses T in feet are determined by equations of the general forn 


y- 
Dy, Fon t Dg ( 
wherein, Cy is a yariable coefficient, which depends upon the discharge rati 
Be ess stream tee" or ai the 4 t tee) and 
Gr 6 p Ga ot Ge ownstream tee) and upon 
whether the flow is divergent & convergent. The velocity head ie be insert 
in Eq. (24) is either om or ee t the upstream tee a or on at the dow 
stream tee) depending upon whether the main conduit flow is directed towa1 
(+) or away (-) from the turbine. In Eq. (16) Cy is replaced by Cp and Ct, f 
divergent flow and by C, and rou for convergent flow. Four sets of coeffici 
are required in order to obtain the piezometric-head changes at each tee fc 
all cases considered. (Table 5 has been utilized for numerical evaluation | 
when D,/D = 2/3.) The coefficients include allowances for changes in veloc: 
head from conduit to riser and vice versa. The details of tee-loss determi 
nation are given in a summarization in Appendix II. . 
With reference to Fig. 1, if the tee losses are ignored, the piezometric | 
heads at points (1), (2), and (7) are assumed to be equal; that is, Py = Po = 
These three heads differ only when the tee losses are introduced. Thus, fo 
example, a 
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| P= PL = 1 wet aks (17) 
Babar se. 3. 2 Oe 
P vy 
=-PevT - she % 

trey bly aaa 2g (18) 


; 

\ both of these equations (applicable to case (a) on Fig. 2) the sequence of 
lriting the subscripts for T is important from the standpoint of avoiding con- 
1sion in signs. Thus, T1,2 indicates Pj - P2 and not Pg - Py. The same 
easoning applies to the downstream tee. 

In all instances of unsteady flow when the tee losses are taken equal to or 
re zero, the total energy heads at the tees are not balanced because of the 
astantaneous difference in velocity heads between main conduit and riser 
between points (1) and (2) on Fig. 2, for example). (The number of times that 
nis would occur during one total solution when the tee losses are evaluated at 
ll other instances during that same solution is relatively small.) This un- 
alanced condition may be assumed to be taken care of by orifice effects in 

1e riser, that is, by the orifice coefficient. If deemed necessary, special pro- 
isions can be incorporated in the machine program to take care of the 
elocity-head differences when T = 0. In this study, such provisions have not 
een made, because over-all results for the Oahe project, at least, would not 
e significantly affected. 


ate-Time and Turbine-Performance Curves 


Transient effects within the hydraulic system are established by turbine 
icket-gate movements. In order for the program to work successfully, the 
ate opening G in per cent must be known at each instant of time during the 
litial interval of governor action, after which the gate opening is assumed to 
emain constant. This statement is only partially true for surge-tank stabili- 
‘studies. (Discussion of what is involved in stability studies is presented in 
nother section of the paper.) The gate-time curve is given by 


G= F(t), (19) 


| which, F,, = a function. It is not necessary that this be a linear function of 
me. Any single-valued function of time may be used, and it may include 
sad and cushion time (for opening and closing of gates, respectively.) Such 
irves are usually furnished by the manufacturer of the turbine. 

In addition, turbine-performance curves are required. They enter the 
rogram as a table involving three variables, whose inter-relationship in 
*neral mathematical form is expressed as 


H = F(G,0,), | (20) 


which F is a function. For carrying out solutions for full-load rejection at 
aximum turbine design head Hmax and for full-load demand at minimum 
rbine design head Hyjn, it will usually be necessary to have tabular values 
G and Qs for values of H well above Hmax and of H well below Hmin. Other- 
ise, the computer will stop because of artificial bounds imposed by lack of 
bular values in regions of computation which the analyst did not foresee or 
‘erlooked. (What this means is that when the computer reaches the end of a 
ble and it has no instructions on what to do next, it stops.) 


46 April, 1959 HY 4 


To perform surge-tank stability studies along conventional lines, the g 
following set of turbine curves is utilized instead of the family of curves indi: 
cated by Eq. (20); 


H= F(a (21) 


5 Bap = Constant 
in which Fq is a function and Bhp = brake horsepower. This equation simply 
gives the relationship between H and turbine discharge Qs5 when the turbine 
power output is constant and can include the effects of varying efficiency. Th 
use that has been made of such curves and the difficulties experienced in the 
computer study are discussed below. 


Machine Computation 


With the aim of keeping machine-programming details to a practical mini- 
mum herein, Fig. 3 has been especially prepared for the purposes of this 
paper. The chart simply depicts the main elements of problem solution but 
not the actual sequence of machine operations.* To show the actual sequence 
of machine operations would require considerable detail, the inclusion of 
which would tend to obscure the main elements of the problem. Fig. 3 applie 
to problem solutions that start from steady-state conditions, but the detailed 
machine program is not so limited. The equation numbers shown in some of 
the boxes refer to the relationships given in Appendix III. It is suggested tha 
the reader first examine that appendix before proceeding with the remainder 
of the paper. 

In box 6 (Fig. 3) all lost heads are computed from the steady-state Q = Qy 
Q3 = Q5 and data given in boxes 2, 3, and 4, and tailwater level is adjusted to 
the value it must have for all heads to balance exactly within the machine for 
the given reservoir elevation ER (box 2). Alternatively, the tailwater level — 
may be kept constant, and ER may then be adjusted for balancing of steady-_ 
state heads. This balancing process insures that the computer starts off witl 
a totally consistent set of values for the variables involved. 3 

Once steady-state conditions are balanced within the machine, a transient 
is set up by applying the initial time interval At (box 5), which fixes the new 
gate setting from G = Fg(t) in box 4. The new G in turn fixes the new H from 
the turbine-performance curves (box 10), using the value of Qs at the begin- 
ning of the time interval. As a result, all other variables are changed, and — 
the machine will compute their values at each step as the transient develops. 
The initial At is arbitrarily selected by the analyst and can be as small as — 
seems warranted. Values as low as 0.01 second have been used to start the 
transient; in most cases, initial At = 0.10 has been used. The computer can 
be instructed through control information (box 5) to use At values of 0.1 up tt 
a total time of t = 1.00 and to use time intervals of one second thereafter up 
to, say, 80 seconds, after which it would use 5-second intervals until the tot 
time for the machine run (given in box 5) is reached. These figures are mer 
ly illustrative. Any reasonable set of time intervals may be used, including 
possible use of a constant At to the end of the run. However, the method of 
solving differential equations (discussed below) will not work expeditiously ; 
the assigned time intervals for step solutions are too large. 4 
*Complete details, including all flow charts and the machine codes, are give 

in the machine program itself, which is being prepared as a separate 
document for distribution on a limited basis. 
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With reference to Fig. 3, computations are carried out in the loop, boxe 
through 25 and back to 7 (but not precisely in the sequence shown), with pr 
out of results in box 26 occurring whenever all the conditions for solution 
each time step are satisfied. The computer may work around this loop mz 
times during any time step before print-out of results occurs. The compu 
may print out all or any number of the quantities indicated in box 26. In bi 
9, the indicated relationship r5 = 0 is, of course, never true. This is sim} 
a trick which is used to save some programming steps for rejection cases 
when G = 0, Qs = 0, and > must equal zero in box 22. It will be noted th: 
the computer determines the head losses (box 13) and tee losses (box 15) ; 
every step. By manual methods, this would entail a tremendous amount 01 
work. To calculate all the hy,’s at each step requires the determination ol 
five Reynolds Numbers and linear interpolations in the tables for six frict 
factors (box 4). Tee losses are evaluated in accordance with all the con- 
ditions specified in Appendix III and by linear interpolation in a table havil 
four sets of coefficients. This degree of refinement in surge-problem co! 
tations is practicable because the machine costs involved are comparative 
small. 

The choice of whether or not to include riser acceleration (box 16) isn 
necessary in future applications. The option was incorporated in this 
program in order to give a direct indication of what effects this item woul 
have on final results. An option on computation of tee losses has also bee 
incorporated in the machine program, but the option is not indicated on Fi 
(box 15). Here again, in future applications, there would be no necessity f 
ignoring tee losses. As long as the form of the tee-loss equations and the 
method of evaluation do not change, any set of numerical values may be ut 
ized for the tee-loss coefficients in the machine program. 


Machine Solution of Differential Equations 


Knowing the values of the primary variables Q1, Q3, Q5, Ws1, and Wsg 
the beginning of a time step, At = tg - t,, the derivatives of these five var 
ables (boxes 14, 22, 23 and 24 in Fig. 3) may be readily determined for ti 
ty. (Qg and Q4 are not designated as primary variables, since their valu 
are easily obtained by algebraic subtractions through use of continuity Ec 
(5) and (6)). The values of the primary variables at the end of the time ste 
tg , are obtained by a process of high-order approximations. (10) ; 

For ease in writing let y equal any one of the five primary variables fc 
the purposes of this section of the paper. Then the values of any one of tk 
quantities at the beginning, mid-point, and end of the time step At are des 
nated as yj, y1,5, and yg, respectively. Roman numerals are utilized to ¢ 
nate approximate values of yy 5 and yg. Thus, y9 indicates the first appr 
mation of yo, yh the second approximation, etc. At time ty, yz and eae 
known. A first approximation is made by utilizing the value of the deriva 
at t, to yield ¥ 


dy. 
oe 1 
Vp = Vy Ot 


Having a first approximation of the five primary variables, as indicated t 
Eq. (22), new values (temporary) of their five derivatives are then compu 


| 
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] 
tly pe Same way as was done at time t1; these new derivatives are desig- 
d dyo . The next step yields a better approximation of yg by making use 


ne ot erage of the two derivatives now at hand, so that 


: i. 
dy dy, 
ual i ab 2 
| af. oath es fae eae de ey 
ling set yu as given by Eq. (23) aside, the computer cuts the original time 
irval in half and determines yj{.5 in exactly the same manner as Yo was 
tulated, with the result that 


; dy 
i ee eee g 
| Y15 = yy+ t pee and (24) 
i igettileay 
II Ae 1 1.5 
eee Ge tea) (25) 


Still another approximation for yg is obtained by applying the same method 
omputation to the second half of the time interval (from t1_5 to tg), so 


aft 
Eitio), 1 15 Ag (26) 
a SE ae ae oes 

II III 

5 a: 
aN time SE igt AS y pee ¥ (27) 
et ar os oA at 


‘final values of all primary variables at time tg are taken as* 


a 
pein daecis 
eee ce 3 (28) 


before these final values are assigned to yg, relative error tests are 
le. The relative error is defined as 


IV ry 
ae Jo Li Yo i (29) 
y Paley: 

3¥5 


| > €,, a pre-assigned margin of relative error, the computer cuts the 
‘inal time interval At in half and goes through the whole procedure (em- 
ied in Eqs. (22) through (29)) all over again. The machine will repeat this 
ess until |Ry| <€,. (For the Oahe project, €, was given a value of 10-3.) 
drevent the computer from taking unnecessarily small time increments, 
machine is instructed to double the original time interval whenever 

3 €9, another pre-assigned margin of relative error. (For the Oahe 

ect €9 was assigned a value of €;/100, that is 10-5.) By the procedure 
outlined, solutions are advanced time step by time step and automatic 

)r control is achieved at each step. 


e factor three in Eq. (28) is based upon theoretical considerations, (10) 
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The necessity for such high-order approximations for solutions of differ: 
ential equations has not been investigated, because the machine time requir 
per total solution is reasonable and because office time could not be spared 
for the mathematical explorations that would be required. It is believed tha 
the machine method of solving differential equations for surge problems, as 
discussed herein, can also be applied profitably to machine solutions ofa 
variety of other engineering problems. 


Application of Machine Methods to Oahe Surge System 


Conditions and Data for Numerical Solutions 


Oahe Dam, located about six miles upstream of Pierre, South Dakota, is 
key project of the system of multiple-purpose reservoirs on the Missouri 
River.(5,6) Each of the scheduled seven generator units at Oahe will have ¢ 
nameplate capacity of 85,000 kilowatts, and hydro-power will be delivered t 
the generators by Francis-type turbines. The hydraulic surge system pert; 
ing to any one of the seven units is schematically shown on Fig. 1. The len 
of the upstream conduit L, is different for each unit due to physical layout. 
Machine results have been obtained only for unit number 6, for which 
L, = 3,820 feet in all cases, except cases Ie and Ie in Tables 6 and 7. The 
two tables and Table 8 give the initial conditions (steady state) and other 
pertinent data upon which manual and all machine computations for full-loa 
rejection, full-load demand, and surge-tank stability analysis are based. 
References to manual computations in all tables and figures herewith apply 
those made independently by engineers of Sverdrup and Parcel, Inc., a firm 
under contract with the Omaha District, Corps of Engineers, for design oft 
Oahe power plant and surge tanks. 

All computations are based upon: surge tank diameters of 70 feet with n 
interconnection between tanks; riser diameters of 16 feet; main conduit and 
penstock diameter of 24 feet; riser height R of 81 feet; orifice diameters of 
12.33 feet; a distance between risers Lg of 40 feet (except for cases Ie and 
Ile); and other data given in Tables 6, 7, and 8. With Lg = 40 feet, the surg 
tanks will, of course, not fit concentrically over the risers. Hence, for cas 
Ie and He the distance Lg was doubled in order to permit a concentric ar- 
rangement of surge tanks and risers, if desired. From a hydraulic stabilit 
standpoint, the distance between tanks should probably be as small as pract 
cable. The eccentric arrangement with Lg = 40 feet has been adopted for tl 
Oahe project. 

All friction factors for both manual and machine computations have beet 
evaluated from data given in Tables 1, 2, and 4. When Np < 105, friction — 
losses for the smallest diameter (16 feet) are negligibly small, and hence, 
values of f are given in the tables for NR< 105, Table 5 has been used to 
evaluate all tee losses for computer solutions. Manual computations in all 
cases are based upon ignoring acceleration effects in the risers and in the 
length L5, but tee losses have been evaluated using Vogel’s data. (4) tty Ii 
noted from Tables 6, 7, and 8 that Ly for manual computations was 3,580 f 
instead of 3,820 feet and that a few other items for these computations di 
slightly from those used for machine solutions. The reason for this is th 
computer runs were made at a much later data, after some Oahe design 


J 
changes had been effected. Other conditions and data as well as results ar 
discussed below. 
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ABBREV |ATED 
TABLE 3 


OAHE POWER PLANT 
TURBINE PERFORMANCE 


Values of Discharge in C.F.S. Versus Gate Opening and Net Head 


“a | a 
LE I* TABLE 2°* TABLE 4+ 6-Gate Net Head on Turbine H in Feet 
a Opening 3: 
. Reynolds Reynolds in & 70 110 | 120 130 | 160 | tas 205 | 240 
Number Number 
& . 2 Ne - 0 0 0 0 0 0 0 0 0 
al R | 3 276 370 | 393 412] 470 | 516 549 | 588 
5 , 0 551 740 | 786 824 | 940} 1031 | 1097 | 1176 
peol ee ses iat ed ROSEN Hee abd 15 sz7 | iio} 1178 | 1237] 140 | 1547 | 16y6 | 1764 
aN aA orale? 108 meine 20 i102 | 1480} 1571 | 1649! 1880 | 2062 | 2194 | 2351 
Here Bal ets igs ae 25 1378 | 1850 | 1964 | 2061| 2350 | 2578 | 273 | 2939 
ae at See an octles 30 1653 | 2220] 2357 | 2473] 2820] 3093 | 3201 | 3527 
ee Aes Renin aie S-O1gs 35 1929 | 2590 | 2749 | 2886| 3290 | 3609 | 3a40| 4II5 
See mae aeoiae mae DiGES 4o 2204 | 2960) 3i42 | 3298| 3760) 4i24 | 4396 | u725 
ea ot me re edt 45 2480 | 3330| 3535 | 3710) 4230 | 4eyo | 4es5 | 5275 
Rees 198 soe 195 DDE? 50 2800 | 3705 | 3910 | 4110] 4665 | 5092 | 5398| 5800 
. : 5 55 3155 | 4065 | 4ze0 | 4y90/ 5090 | 5532 | sa¥0| 6260 
1 60 3460 | 4430 | 4655 | 4870| syes | so¥s | 6278] 6740 
65 3745 | 4775 | solo | 5230| says | 6322 | 6esz | 7235 
70 oso | 5135 | 5375 | s600| 6215] 6692 | 7050) 7585 
Von Karman-Prandt] Curve for hydrodynamically smooth 75 4355 5430 5670 5900 6545. 7048 7428 | 8015 
aries. 80 4600 | 5750| 5985 | 6215! 6875 | 7390 | 7778| #8330 
85 4900 | 6035 | 6285 | 6525) 7190| 7705 | 8095) e715 
t 90 5195 | 6280| 6530 | 6770 7475| s025 | ay25| 9015 
d upon Relative Roughness e/D = 1.25 x 107 95 5450 6490 | 6735 6975| 7685 | 8250 8677 | 9355 
: 100 5705 | 6670| 6910] 7160/ 7875| e458 | 8905] 9665 
id upon Relative Roughness e/D = 3.3 x 107 
ee 


es in Tables 2 and 4 taken from U. S. Bureau of Reclamation 
graph #7, Figure 4, 


Adapted from Omaha District Curves, Mode] 734-H, dtd 20 November 1957. 


Discharge Ratio 


conduit 


0.2 -0.47 

0.3 -0.60 

0.4 -0.69 

0.5 -0.75 

0.6 -0.75 

Z. 0.7 -0.77 
‘ 0.8 -0.75 


— 


TABLE 5 
Tee Loss Coefficients” 


0, /o = 2/3 


Diverging Flow 


Converging Flow 


0.92 0.68 0.83 0.15 
1.50 0.86 1.12 0.26 
2.16 0.97 1.38 0.41 
2.4 1.07 1.62 0.55 
3.75 1.10 1.82 0.72 
4.77 1.10 2.00 0.90 
5.95 1.05 2.17 1.02 


+ 


Use For 
Tee Losses 


Ty ,2 T,9]71,778,14 


* Prepared by Omaha District, 
** Discharge ratios may be either Q7 , 


CE, from Vogel's data (4). 
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its for Full-Load Rejection 


{th reference to Table 6, all cases are based upon rejecting full load on 

irbine, 142,000 horsepower output, at a maximum net head on the turbine 
feet for steady state. For this condition the gate opening is restricted 

}3 per cent for all computer solutions (70 per cent for manual compu- 

4s) in order to avoid overload on the generator, and closure of the gates 

hieved in 3.31 seconds (computer) at a uniform rate of 20 per cent per 

ad. Turbine-performance values (Eq. (20)) have been obtained by linear 
lation in a table of the form indicated by Table 3. The actual turbine- 

jrmance table that was used in machine computations gave Qs, values for 

i H at closer intervals (AG = 5 between 0 and 100, AH = 5 feet between 

id 240 feet). 

gs. 5 and 6 indicate some of the results that were obtained for the re- 

om cases. Fig. 5 depicts the water surface elevation of the downstream 

2 tank and the difference in water surface elevation between the two tanks 

function of time for all cases indicated in Table 6. Case I and Manual 

outations are based upon the same data and conditions, except that for the 

ial computations Ly is over 200 feet shorter, initial Q5 and G are slight- 

ther, and acceleration heads in the risers and in the length Ls are not 

uted. The basis of computation for all other cases shown on Figs. 5 and 

he same as for Case I, except that for: 


ase Ib, acceleration heads in the risers have not been computed. 

ase Ic, the orifice coefficient is 0.8 instead of 0.7. 

ase Id, the orifice coefficient is 0.6 instead of 0.7. 

use Ie, the length Lg is 80 feet instead of 40 feet, Ly and L5 being some- 
what smaller as a consequence. 


‘amination of the curves on Fig. 5 indicates successful application of 

ine methods to rejection problems. It is noted that: (a) the numerical 
of the orifice coefficient has the most pronounced effect upon the water 
ce elevations in the tanks and (b) in comparison with the results of manu- 
mputations, the difference in water surface elevations between tanks is 
antially greater for all computer cases during the first quarter cycle. 

1 machine-computation time for each of the cases shown on Fig. 5 aver- 
slightly over one hour. The Q values for Case Id on Fig. 6 are more or 
typical of those for all other computer cases indicated on Fig. 5. 


its for Full-Load Demand 


ith reference to Table 7, all cases are based upon picking up load when 
sservoir elevation is at minimum design level (ER = 1540). The rate of 
ng the wicket gates is linear in all cases, the rate being 20 per cent per 
d (governor time of 5 seconds) for all except Ila. For Case Ila, the 

nor time is 10 seconds. Applicability of Table 3 to load-demand cases 
-same as for the load-rejection cases previously discussed. Case II and 
al Computations (Table 7) are based upon the same data and conditions 

t that for the latter: L 1 is over 200 feet shorter; initial steady-state G 
-o (instead of 5.0); gate opening at end of governor action is 90 per cent 
ad of 100 per cent); tailwater elevation is 1420 (instead of 1421.7); and 
eration heads in the risers and in the length Ls are not computed. All 
iter cases are based upon an initial steady-state discharge Qs of 390 
(corresponding to initial G = 5.0) instead of zero. A comparatively small 
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finite value of Q5 at the start seems to be necessary because of the machine 
methods of solving differential equations that has been utilized. A smaller 
value than 390 c.f.s. might have been used successfully. However, provided 
that the initial Q, is small in comparison to the full-load Qs, the over-all re: 
sults are not likely to be significantly affected. 

Figs. 7, 8, and 9 indicate some of the results that were obtained for the 
demand cases. All data shown under Case II in Table 7 also apply to Cases 
II-1 and II-2, except that the surge tank base elevation E, of 1515 does not aj 
ply to II-2. Under Case II-1, the surge tank water surface drops down into 
the riser, but under Case II-2, Eg is sufficiently lower (about 3 feet) so that 
the water surface does not recede into the riser. (With this small difference 
the change in R is ignored.) The basis of computation for all other computer 
cases is the same as for Case II-2 except that for: 


Case Ila, the governor time is 10 seconds instead of 5 seconds, and 
with E, = 1515 the water surface drops down into the riser. 

Case IIb, acceleration heads in the risers have not been computed. 

Case Ild, the orifice coefficient is 0.6 instead of 0.7. 

Case Ile, Lg is 80 feet instead of 40 feet thus causing Ly and Ls to 
be somewhat smaller, and K, is 1.00 (re-entry tube) instead of 0.5. 


Comparison of manual and machine results on Fig. 7 demonstrates that 
machine methods have been applied successfully also to full-load demand 
problems. The question of whether or not to permit the water surface to dre 
‘down into the riser in actual operation is immaterial insofar as the purpose 
of this paper is concerned. Results are now available as a basis for decisioi 
in that respect. It will be noted from Figs. 7 and 8 that machine results indi 
cate oscillations of rather large amplitude between tanks for a short period, 
when the water surface is permitted to recede into the risers. Cases II-1 at 
Ila are directly comparable except for governor time, and hence results for 
these two cases are plotted separately for comparison on Figs. 8 and 9. : 
Doubling the governor time (IIa) seems to make little difference in over-all 
results for the Oahe surge tanks as such; that is, low point of water surface 
elevation in the risers (Fig. 8) is about the same and the timing therefor is 
also about the same. But, of course, with a slower governor, pickup of the — 
turbine load is materially delayed, as comparison of Qs values on Fig. 9 ind 
cates. It is interesting to note from this figure that except for a slight lag 
time, the results of Case Ila are almost a carbon copy of those of Case II-1. 
This demonstrates that the computer does the job with slavish adherence to 
instructions, and, as a consequence, the analyst cannot afford to slight or 
overlook any pertinent detail in the preparation of the machine program. A¢ 
al machine-computation time for each of the demand cases discussed above 
averaged about one and one-half hours. 


Surge-Tank Stability Analyses 


In normal practice, stability analyses are made for small power-demand 
changes at minimum design head on the turbine in the region of drooping ef- 
ficiency (between gate opening for maximum efficiency and full gate.) Ideall 
what is desired is to go directly from the initial steady-state power Bhp; to 
a new steady-state power Bhpj. For example, on Fig. 11 it is desired to pr 
ceed directly from the point Mj (Bhp; = 75,000, Hy = 120) to point Mo 
(Bhpo = 78,000, Ho = 119,6), the new steady state for constant turbine output 
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of 78,000 horsepower. However, this is not possible, for sudden gate 
movements in response to rapid governor action set up transients in the w 
surge system. These transients will die out ina reasonable length of time 
the system is properly designed. By conventional methods, it is assumed | 
purposes of analysis that once entry onto the curve Bhpg is made, at whate 
point that may be, the power delivered by the turbine will have a constant 
value Bhpg thereafter. The oscillations in Qs, H, and G which occur with 
Bhpg constant are then determined by manual methods of arithmetic inte~ 
gration. (2) Such methods are necessarily based upon simplifying assumpti 
as a result of which the total transient problem is treated as one of separa 
multiple parts. In essence, the hydraulic system is divorced from the inte 
connected mechanical and electrical system for the practical purpose of ai 
lyzing the hydraulic transients as such. Use of this same approach to so- 
lutions by computer (employing the method of solving differential equation: 
previously discussed) has been unsuccessful, but only in the sense that the 
computer has failed to produce the same results as those obtained by mani 
computations. In another sense the computer study has been quite fruitful 
that it has forced a reappraisal of the total “transient” problem and has hi 
lighted the shortcomings of manual methods of analysis. Discussion of the 
need for comprehensive analysis of power plant transients as an integrabl 
whole is presented in another section of the paper. 

With reference to Table 8, all computer runs were made starting with 
Bhp, = 75,000 for initial steady state at a net head on the turbine of either 
or 130 feet (the former value being the minimum design head for the Oahe 
turbine.) The gate-time curves indicated in Table 8 are shown on Figs. 4 
10, on which Go is the initial steady-state gate opening (line one of Table { 
Curve No. 1 on Fig. 4 represents part of the full-range gate-time curve 
furnished by the Oahe turbine manufacturer. Its relation to Curves No, 2, 
and 4 is indicated on Fig. 10. Curves 2, 3, and 4 are entirely arbitrary. | 
for all computer runs are exactly the same in all respects except for indi 
ed initial H values (line 2 in Table 8), gate curves, and indicated inclusion 
exclusion of acceleration heads in the risers (line 20 in Table 8). 

The first computer run that was made for stability analysis of surge ta 
has been designated Test No. 1 and involved use of Gate-Time Curve No. 
(Table 8). Data for this test are exactly the same’as for Manual Compu- 
tations, except for the inclusion of acceleration heads in Ls and the differ 
in Ly. The path taken by the computer in trying to reach Bhpg = 78,000 is 
indicated by Test No. 1 on Fig. 11. Full gate was reached at 1.1 seconds’ 
fore Bhpy was reached. In following the paths indicated by Tests No. 1,2 
and 3 on Fig. 11, the computer was instructed to follow the gate-time cur 
given to it (in tabular form) until it reached the curve Bhpg, after which c 
trol for development of the transient would revert to that curve. (That is, 
thereafter the function represented by Eq. (21) instead of Eq. (20) would h 
to be satisfied as well as other relationships.) Also, the computer was pI 
grammed to stop whenever it reached full gate, since oscillations past ful 
gate are unacceptable. 

Test No. 1 having failed, Test No. 2 was started from a region (H = 13 
where the efficiency curve is flat and turbine efficiency is practical const 
down to H = 120. Again, the computer followed a devious path, reaching E 
at 1.36 seconds (Fig. 11), after which it continued along the curve Bhpg tc 
gate. With the thought that, perhaps, riser acceleration effects caused th 
devious path for Test No. 2, Test No. 3 was started from the same point 
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1130) with the effects of riser acceleration excluded. This time, the path 
ed by the computer was satisfactory until the Bhpg curve was reached 
5 seconds (Fig. 11), but after that the computer followed the Bhp»s curve 
ird to full gate once again. These failures to establish a transient (simi- 
p that obtained by Manual Computations) that would produce oscillations 
5, G, and H along the constant Bhp curve of 78,000 without at any time 
in full gate do not indicate that the Oahe surge system is unstable in 
egion investigated. Rather, they indicate that any transient set up in the 
er of conventional methods of analysis is likely to resemble the actual 
et only by accident, for these methods ignore (with justification*) the 
lations of the gates with time due to other effects at the turbine and 
mnor end of the surge system. In anticipation of any doubts that may arise 
erning the adequacy of the machine program, sufficient data are given in 
paper for the reader to assure himself that the computer has done the job 
}cordance with instructions. Machine print-out of results for Test No. 2 
tiven in the upper part of Table 9. These data together with the equations 
ppendix III, data given in Tables 1, 5, and 8, and gate-time Curve No. 1 on 
4 will permit checks on machine arithmetic in a relatively short time. 
ll be noted that the machine does not print out values of hao and hg 4B they 
be obtained by separate computation. 
hen a load change is called for by external electrical power require- 
‘S, oscillations are set up in the turbine governor system itself during the 
several seconds. By virtue of interconnection, the gates also oscillate. 
. effort to represent at least the first one-half cycle of these oscillations, 
-time Curves No. 2, 3, and 4 (Fig. 10) were invented, using the “dead” 
“cushion” ends of the manufacturer’s curve (of which Curve No. 1 isa 
to form the peak regions and the tail-ends. Application of these gate- 
curves then gave the results shown on Fig. 12. In all three cases the 
juter started and stayed within a region in which turbine efficiency is 
tically constant. In following the path indicated by Curve No. 2 on Fig. 12, 
omputer reached Bhpg at 0.8 second, after which it proceeded along Bhpg 
ll gate in the same manner as before. (Machine print-out of some of the 
its based on use of Curve No. 2 is given in the lower half of Table 9). At 
point in the study a program change was made by which the machine was 
ucted to ignore the first intersection with Bhpg and to follow gate-time 
e No. 2 (Fig. 10) until it reached Bhp a second time, after which control 
d revert to Bhpj. For this condition the machine stopped at 1.2 seconds 
‘the end of gate-time Curve No. 2 shown on Fig. 10), and this test was 
joned. This same procedure was applied to gate-time Curves No. 3 and 
h the results shown on Fig. 12. In following the path marked Curve No. 3, 
omputer followed instructions perfectly, reaching Bhpg a second time at 
} second, after which it followed Bhpg downward to the right and stopped 
[25 seconds. It stopped at this point because it had reached the end of 
ible giving H as a function of Qs along the curve Bhpg, the end of the table 
g been purposely set at a limit of H = 140. Gate-time Curve No. 4 was 
ed with the aim of having the path loop back and intersect Bhp» in the 
/ 
ification exists because, normally, these oscillations cannot be described 
juately without bringing into the problem all the relationships that tie the 
-aulic, mechanical and electrical systems together. Manual analyses on 
basis would be totally impracticable. 
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opposite direction from that which had been indicated by Curve No. 3. The 
computer started a return loop along Curve No. 4 as indicated by the points 
at 0.9 and 1.0 second on Fig. 12 and then stopped at 1.0 second (the end of 
gate-time Curve No. 4 on Fig. 10). At this point it was decided that there w 
altogether too much arbitrary judgment being exercised, and further compu 
er studies of stability were not undertaken. 

Machine-computation time for all these tests was comparatively small. 
However, a good deal of time was spent in analyzing machine results and in 
trying to establish a transient using a function of the form given by Eq. (21) 
as acontrol. It is concluded from the experience thus gained that analyses 
utilizing this procedure yield artificial results and that, as a consequence, 
the results of conventional methods of analysis of surge-tank stability along 
the lines previously described are incomplete and inconclusive. This does 
not raise or imply any reflection on the excellent work that has been done i 
connection with hydraulic design of power plants. For, it has long been rec 
nized that simplifications are a practical necessity in order to render powe 
plant “transient” problems tractable for solution by manual methods in a 
reasonable length of time. However, successful adaptation for computer so 
lutions of very complex problems in many fields of engineering and science 
now raises the question as to the necessity for continuing to use simplified 
procedures that yield results that are not entirely satisfactory. 


Some Ideas for Comprehensive Study by Digital Computer 


The ideas presented in this section of the paper have been extracted frot 
an unpublished “Proposal for Comprehensive Study by Digital Computer of 
Hydraulic and Governing Transients Affecting Design and Operation of a 
Power Plant System”, prepared by the writer on his own initiative and sub- 
mitted to the Corps of Engineers for consideration. This proposal includes 
the derivation of necessary equations and suggested methods of machine ani 
sis and computation for undertaking a comprehensive study, which is only — 
briefly discussed below because of space limitations. The only reason for 
viting attention to the proposed study at this time is to indicate to the engi- 
neering profession that efforts are being made to overcome the limitations 
conventional methods. ‘ 

Formidable limitations of time and means have required up to now that ; 
power plant system be treated as a divisible whole for purposes of analysis 
of hydraulic and governing transients induced by external power-load vari- 
ations. (As used herein, the term “power plant system” includes the surge 
system as previously defined and one power unit, turbine and generator, su 
plying load on an independent basis.) With the advent of powerful, high-spe 
computers, it appears that, for the first time, it is practicable to analyze tl 
total “transient” problem as a “feedback” system of interwoven effects ratl 
than as a problem having separable, multiple parts. Studies of water hamn 
surge-tank effects, turbine-speed regulation, and governing stability need r 
longer be made on an essentially separate basis if a computer program car 
be developed to handle the total problem as such. Recent analog-computer 
studies by Koenig and Knudtson(9) are only a step toward the ultimate goal, 
for these studies involve simplifications which tend to cast some doubt on t 
validity of results obtained thereby. d 


SURGE TANKS 67 


| appears that studies of the “total problem” can be made with conclusive 

re only when the relationships controlling hydraulic behavior in the 

uit, penstock, and surge-tank system are solved simultaneously with the 

‘ecting relationships between turbine output, the governor, machine in- 

a and speed, and external power load. To do so requires the simultaneous 
ion of a great number of equations, short time step by short time step. 
time step would probably have to be in the order of hundredths of a 
nd (or less in some cases) in order to account for all pressure dis- 

ces and reflections* generated during the existence of a transient. 
ously, such an approach is totally impractical, if not impossible, for 

ual computation, but, apparently, not for a modern digital computer hav- 

the capability and computational speed comparable with those of an IBM 

computer. 

. has been roughly estimated by the writer that the machine time required 

total solution on such a computer would be well within practical range, if 

1ecessary machine program can be successfully developed. A total so- 

yn would yield time histories of velocity and piezometric heads at a large 

ber of pre-determined points within the surge system as well as time 

ories of turbine power output, turbine discharge and head, turbine and 

‘rator speed for any assumed or given external power-load variation with 
‘ applied at the power unit under study. Important quantities such as 

fe tank size, orifice characteristics, penstock and conduit characteristics, 

hine inertia, governor characteristics, etc., would all be introduced in 

otal problem as parameters. By the proposed methods, water-hammer 
mass-oscillation phenomena would not be treated as separate entities. 

.e from the basic assumption of essentially independent operation of the 

sr unit, no other simplifications of important significance appear to be 

‘Ssary for the methods of solution that have been proposed. To accomplish 

his on a digital computer requires the simultaneous solution by numerical 

10ds of partial-differential equations within the surge system and of ordi- 
differential equations at the turbine end of the system. Methods for do- 
his are apparently available. It remains to be determined whether or not 
undamental partial-differential equations describing hydraulic behavior 
in the surge system can be solved by computer on a practical basis. The 
er believes that this can be done, but a machine programming job of great 

e will have to be performed before success is achieved. 


CONCLUSIONS 


machine program for the IBM 650 computer is now available for solving 
-rejection and full-load demand problems of a surge system comparable 
at at the Oahe project in one to two hours of machine-computation time 
“ase. The machine time per solution for installations having only one 

- connected to the main conduit will be materially less. The program de- 
ed herein may be easily adapted for such cases. The relatively small 
required per machine solution and the comparatively small costs in- 

1d make it feasible to enlarge the scope of design studies of surge 

ms well beyond the limits ordinarily imposed by manual methods of so- 
n. By examining the results of a number of machine solutions readily 


1e study that is envisioned, effects of the compressibility of water and the. 
ticity of flow boundaries would be included. 


, 
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obtained as a result of changing one or two major factors from one machine 
run to another, a better understanding of problem phenomena for a particula 
installation will be obtained. As a consequence, construction savings may bt 
achieved in some cases. 

A digital computer program which will yield solutions for transient effect 
resulting from small-load changes from one steady-state power output to an 
other is not yet available despite considerable efforts to develop such a 
program. It may be entirely possible to develop such a program by utilizing 
methods of solution other than those described herein. However, it appears 
that any method which ignores the applicable relationships between turbine, 
governor, machine speed and inertia, and external power load for the purpos 
of analyzing the hydraulic transients as such is bound to yield inconclusive 
results. The premise that such simplifications give conservative results in 
sofar as surge-tank stability is concerned may be true but the degree of con 
servatism is not known in application to particular cases. For a long time < 
need has existed for analyzing the total “transient” problem without the maj 
simplifications that are ordinarily made. When only manual methods of ana 
sis and computation were available, fulfillment of this need on a practical 
basis was not possible. With the advent of powerful, high-speed computers 
the prospects are now excellent that this need can be wholly satisfied. Some 
efforts have already been made toward the ultimate goal of adapting the tota 
“transient” problem of a power plant for solution bya digital computer in su 
a way that machine results for water hammer, surge-tank effects, turbine- 
speed regulation, and governing stability would be delivered in one package | 
a few hours of machine-computation time. If successful adaptation is ulti- 
mately achieved, many benefits and possible construction savings will accru 
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APPENDIX I 

Notation q 

A - Cross-sectional area in sq. ft. of conduit and penstock. 
Ay - Area of orifice in sq. ft. : 


A, - Area of riser in sq. ft. ee 
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- Area of surge tank in sq. ft. 


) 
: 
ap - Brake horsepower (output of turbine) 


L - Coefficient of orifice in Q = CgAg V2gho 
: t 

: - Coefficient for risers, a ree 

i 2gar 2gA2 


t 
D» Cp - Coefficients of loss for diverging flow in tees, function of dis- 
| charge ratio (Qyiser/Qconduit) 


hs C,. - Coefficients of loss for converging flow in tees, function of dis- 
charge ratio (Qriser/Qconduit) 9 
- General coefficient of loss for the tees in T = Cy oe 


| - Diameter in feet of conduit and penstock. 
- - Diameter of riser in feet. 
; - Diameter of orifice in feet. 
- Diameter of surge tank in feet. 
- Absolute roughness length in feet for concrete. 
- Absolute roughness length in feet for steel. 
- Pre-assigned margin of relative error. 
2 - Elevation of reservoir in feet above M.S.L. 
- Elevation of base of surge tanks in feet above M.S.L. 
- Elevation of centerline of turbine distributor. 
- Darcy-Weisbach friction factor, function of Reynolds number. 
- Friction factors in upstream conduit. 
,f{, + - Friction factors in upstream and downstream risers, respective- 
ly. 
- Friction factor in penstock between risers (Length Lg). 
- Friction factor in penstock (Length Ls). 
. Fg - Functions of Reynolds Number for concrete and steel. 
Fg, F,- Functions defined by Eq. (111-46), (111-45), and (I-49). 
- Acceleration due to gravity in ft. /sec.2 
- Gate opening in per cent for turbine wicket gates. 
- Initial gate opening for steady state. 
- Surface resistance loss in feet. 
- Form loss in feet. 


- Head loss in feet (may include all types of losses and velocity 
head). 
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- Head loss in feet in upstream conduit in length Ly (includes 
all form and surface-resistance losses as well as velocity 
head), defined by Eq. (I-21). 


hyo» hy, - Head losses in upstream and downstream risers, respective. 
defined by Eqs. (III-26) and (III-29). 


hy, 


hy, - Head loss in penstock Lg, defined by Eq. (I-24). 
hy5 - Head loss in penstock length Ls, defined by Eq. (II-25). 
— Net head on turbine in feet. 
ho - Head on orifice in feet (piezometric drop). 
hg - Acceleration head in feet, defined as i 2 
hay - Acceleration head in upstream conduit (Lj). 
hao Rag - Acceleration heads in upstream and downstream risers, re- 
spectively. 
has - Acceleration head in penstock length Lg. 
a5 - Acceleration head in penstock length Ls. . 
hy - Velocity head in feet. 
K - A constant, in general. 
Ky - A constant defined by Eq. (III-22). 
K} - A constant, $e x 
Ko, K4 - Constants which apply to upstream and downstream risers, | 
spectively, defined by Eqs. (III-28) and (I-31). 
Ky - Entrance loss coefficient for upstream conduit. 
Ky - Entrance loss coefficient, surge tank to riser. 
Ki; - Loss coefficient for trashracks. 
K, - Loss coefficient for all gate slots. 
Kts - Loss coefficient for all transitions. 
Khb - Loss coefficient for all horizontal bends. | 
Kyp - Loss coefficient for all vertical bends. : 
Kot - Coefficient to account for all other form losses. q 
Ky - Velocity-head coefficient = 1.0. : 
Lj, L3,L5- Lengths in feet in main conduit (see Fig. 1). 
Le - Portion of Ly lined with concrete. : 
Lg - Portion of L, lined with steel. 
Np - Reynolds Number, ¥D 3 
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P - Piezometric head in feet above M.S.L. 

Py, Ps, PP... - Piezometric heads at indicated points (see Fig. 1). 

Qi - Discharge in c.f.s. in upstream conduit. 

Qs, Q4 - Discharges in c.f.s. in upstream and downstream 

) risers, respectively. 

Qs - Discharge in c.f.s. in penstock length Lg (Fig. 1). 

Qs - Discharge in c.f.s. in penstock length Ls (Fig. 1). 

le| - Indicates absolute value of Q (ignore sign). 

R - Height of riser in feet (Fig. 1). 

Ry - Relative error of y variables. 

11, 13,15 _ - Constants defined by Eqs. (II-1), (III-2), and (II-3). 

ro, r4 - Constants defined by Eqs. (III-4) and (III-5). 

t - Time in seconds. 

At - Time interval. 

po - Tee loss in feet. 

Ti2 - Tee loss in feet between points 1 and 2 (Fig. 1). 

T1,7> Tg,9, T8,14. - Tee losses in feet between indicated points. 

Tw - Tailwater elevation in feet above M.S.L. 

Vi - Velocity in ft./sec. in length Ly. 

Vo, V4 - Velocities in ft./sec. in upstream and downstream 
risers, respectively. 

V3 - Velocity in ft./sec. in length L3. 

V5 - Velocity in ft./sec. in length Ls. 

y - Used to represent any one of variables, Q;, Q3, Qs, 
Ws, and Wsp. 

Ws} - Water surface elevation (M.S.L.) of upstream surge 
tank. 

NS5 - Water surface elevation (M.S.L.) of downstream surge 
tank. 

y - Kinematic viscosity of water in feet squared per 
second. 

Sp . - Functions defined by Eqs. (I-15) and (I-16). 


,W - Functions defined by Eqs. (III-18) and (I-19). 
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APPENDIX III 


Mathematical Relationships 


Notation 


Detailed notation is given in Appendix I. In general, the following notation 
is utilized, based upon the pound-foot-second system: 


A 
Cc 
D 
E 
e 
F 
f 

G 


- Cross-sectional area H_ - Net head on turbine 

- Coefficient (riser, orifice, tee)K - Coefficient for losses 

- Diameter L_- Length in conduit 

- Elevation (M.S.L.) NR - Reynolds Number 

- Absolute roughness length P_ - Piezometric head (M.S.L.) 
- Function Q_ - Discharge 

- Friction factor R_ - Height of riser 

- Gate opening in per cent T - Tee loss 


. 
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a - Acceleration head t - Time 

jz, ~ Lost head V_ - Velocity 

e - Velocity head Ws - Water surface elevation 
: 

} 

. Convention 


“he net head on the turbine and all piezometric heads are always positive. 
signs of all lost heads and acceleration heads are automatically assigned 
e€ equations that define them. The signs of all velocities will be the same 
ose of the corresponding discharges, since Q= AV. The sign convention 
all discharges is as follows: 


All flows Q directed toward the turbine are given a positive (+) sign. 
All flows Q directed away from the turbine are given a negative (-) sign. 
_ Flows Q in the risers are positive when directed downward and negative 
when directed upward into the surge tanks. 


hematical Relationships 


Che physical equations and functions required to solve surge problems for 
system shown in Fig. 1 are given in this appendix. The major variables 
ch vary with time t are: all Q’s, all V’s, all P’s, all -f’s, Ws,, Wsg, G, H, 
Cp, Ce, and Ce Such quantities as As, R, Eg, D, Co, etc., may enter the 
shine program as parameters; that is, such quantities are constant during 
One run, but may be varied from one run to another. 


ag | z 
<p so Abe eae Midis og ore (111-1) 
dQ3 as poe Palen r epee i) 
ae er3( a 3 L,) 3 13 (11-2) 
dQ5 _ e eee 
Her = BF5(Py, - Pus - by) eo (II-3) 

A 
6 To += , Ws, 2 E, 
E = &tp(Ws, - Py - hy ) € my) 
2 
: Ay 


Poy ea ee 
2 We, eR? 1s Bs 


Ar 
Tha pe We Sue. 
: = gr), (Ws, - Py - a << . (II-5) 


ry fee ee ee ee 
4 Wegoeiaon 222 Sha 


%+%=% (III-6) 
B#QaS > oF YHG+Q—o we 
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aay | dae _ 403 


+ 2, . eo ~ Le 
Sasa tye (by differentiation of III-6) (1-6 
dQ3 dQ  dQ5 dQj dg@ dQ das 
ata aerate ae at (m-$ 


Substituting from (II-1), (II-4), and (IM-2), Eq. (III-8) becomes 
ry(E, - Py - by ) + 7p (Ws) - Pp - by ) = 7, (P, - Pg - hy ) (111-10 
1 2 3 
Substituting from (III-1), (III-4), II-5), and (II-3), Eq. (III-9) becomes 
hes eer + Pe Sal ag +e (iparFomhe } agg ay (I-11 


Pie he ea a ee ees (I1t-12 


Pg a Po=18 9 and Pa = Ply Tg ah (I-13 


Using (III-12), eliminate Py and Py from (I-10) to obtain 


Pg= a P| +f, in which an 


ry +o +73 a 
ie r3 —, an (mmt-1! 
ee iy -h a Sta tana mrs) 6 
‘a 1,7 “Ly ¥, RO Ly FO 1 “1,2 in" (miI-16 


Using (III-12) and (I-13), eliminate Py, Pg, Py4 from (I-11) to obtain 


Pa = re) Ft y, in which (11 
-r]1 - r2 ) 
= 7, + 5 » and (m-16 

fag a peek en 
y= T+ Fs (z, = ty +a Ts (Ws)+T, 5 Ly) | 
(I-15 

rh Ts 
t+ fs rP = (Ws, + Tg,9 - by.) + + 55 ze 5 (Py5 a Tg, a) 
Solving (II-14) and (III-17) simultaneously yields 

y- 4 
P, = deed (1-26 


a-O6 

Knowing P, from (I-20), Pg may be obtained from either (I-14) or 17. 
Pg, P7, Pg, and Py4 may then be evaluated by use of (III-12) and (IMI-13). Te 
evaluate (III-20), it is first necessary to determine the head losses hy, and te 
losses T in (III-16) and (I-19). ; 


RE goer } 

Head Losses are determined by Eqs. (III-21) through (III-31), as follows: 
m= [9 Gtr AE]] [a “(2 
K = See with So ETE Cae aE ig SERIO 18 SB Te Me (m2 
Ks : 2 - 
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e 


c 
= constant 
e 


f£, = FN) =F GP), 


VD 
) fos P(N.) = EAs) ; = = constant 
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(I-23) 


Pables of values for f, and f, are needed for machine runs; see, for ex- 


ple, Tables 1, 2, and 4 herewith. Similarly, for all other friction factors, 


ests aS 


my = Kf, 2, |e.) a, 


) 
| A CR ie 
i L o Ke 2 Dee | 25] & 
2 2gc A 
i ome) 
: al 
a es") when Q is negative and Ws, =5.. 
eA. 
ie K, 
oe = ae when Q, is positive and Ws) > «&E.. 
2A 
Ey 
a = O, when Q, is positive and Ws, = Eo 
C.. = 0, when Ws) SES, irrespective of sign of Q,° 
K, = a 5 » when Ws, = E- 
2g) A. 
Ws, -E +R 
K, = ——jF—, when Ws, < E 5: 
2eD_ A. 
a (c, +2, +s) [%] % 
z 2gCKA\ 
[oMie) 
Cc = when is negative and Ws eS Ee 
ree 2? qi, 2S 
2gA 
aB ‘ 
t+ K, 
Cas 5 » when Q, is positive and Ws, > E.- 
2gA. 
C.. =) 0, when , is positive and Ws, = Es: 
CL. = 0, when Ws, < Eo» irrespective of sign of Q 
K, = x 5 » when Ws, 2 E.. 
Soe 
Ws, - E_+R 
Ky, = = » when Ws, <4 Eo: 
2gD A 


Fg, £4, f5, all of which are functions of NR, and same tables may be used. 


(111-24) 
(I-25) 


(IIT-26) 


(I-27) 


(I-28) 


(111-29) 


(III-30) 


(10-31) 
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Tee Losses are evaluated in accordance with the following summarization 
Cases A through E (Eqs. (III-32) through (III-39)). The coefficients Cp, Cp, 
Ce, and Ca are evaluated from tables; see, for example, Table 5 herewith, 


For ease in writing, let hy; = ee hy3 = oe and hys5 = os and recall tha 
Ty,2 = Py - Pa, T1,7 = Pi - Pv, etc., as indicated by Eqs. (III-12) and (II-1i 
Case A - Tee losses ignored. 
P, = Pg = P7 and Pg = Pg = P14 


Case B - Q, and Q3 are both positive and Qg is negative. 
Q3 and Qs5 are both positive and Q,4 is negative. 


52 =; Cy hy Tg,9 a Cp me (1-3! 
le iS) 
T ~ =C'h Te 47> Cy h (111-3: 
ley) Vy 8,14 D we ; 
Q 

Cy and C, = function a Cy and Cy = function a Table 
Case C - Q, and Qg are negative and Qg is negative. 
Q3 and Qs are negative and Q, is negative. 

To = (Cy - Sp) "ys Ty 9 = (Cy - &) ™s (1-3 

77 es: Ty, 1h =p *V. (1-3 

9 : 

Cy and Ch = function = C, and Ch = function a Table 

Case D - Q, and Q3 are positive and Qg is positive. 

Qg and Qs are positive and Q, is positive. ¢ 
7 CMa at ly Pee. ' 

1,7  “e v3 8,14 = Cc) hee (m1-3 

ie r (cy ff ae Tg,9 cs (c, 3 =) AP (1-3 

C, and Ch = function a C, and Cc! = function be Table 

3 ° ee | 
Case E - Q; and Q3 are negative and Qo is positive. 
Q3 and Q5 are negative and Q, is positive. 

fT. 5 = <0h Ty ig #-09 

Lae C0 8,9 eC V3 (11-3 

7, _ =-C! : 

fy ax OG le intr (mI-3 

(04 Vv. ' 

% oie 

C, and C! = function a. C, and C! = function et Table 


As long as the form of the equations and the method of evaluation of tee 


losses do not change, any set of numerical values may be used for the 
coefficients in the machine program. ! 


a 

i 
+ Callie a, eA 
| RaQ, 
fe cA at 

; 

ae 

gA 


tn r a 


si 
a 
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; h = 
a 
5 


gA 


cceleration heads are defined as follows: 


as: 


at 


R, when Ws, 2 E. 
at s 


Ws) - E. + R, when Ws 


1 


R, when Ws, => E, S 


Ws, - E. + R, when Ws, < Eo 


<E. 
s 
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(111-40) 


(I-41) 


(111-42) 


Water surfaces in surge tanks depend upon following equations as well as 


iers: 

| Ss 

& A 
ty 

t el a 
ty 


Aes when Ws) > E.: 


Aas 


when 


Ws 


2 


= E, and Q,, is 
+5. and Q, is 
< E.- 
> E.- 
= and Q), is 


= E. and Q, is 


< Eo: 


" (I-43) 


- (I-44) 


At Turbine end of the surge system; the following relationships must also 


satisfied: 


G = F(t) 


(11-45) 


sives the gate opening as a function of time. Either given by manufacturer 
nust be assumed. The function enters the machine program as a table; 


ce, it may have any form. 


H = F(G, Q ) 


(11-46) 


rurbine-performance curves enter the machine program as a table; see 
example, Table 3 herewith. In program, G and Qs are the independent 
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* 


variables and H is the dependent variable; that is, H is computed when G 2 ‘ 


Qs5 are known. : < 
oe) 
ey (11-47 
a= Pis + Sake Tw, in which ‘ T-4' 
Mae alae) fearer dee. (I-41 
é' a 2 
H (In1-4! 


Fy (@pnp = constant 


This function used only for stability studies; see, for example, Figs. 11 2 
12 herewith. These curves enter the program as tables. 
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_ TWO METHODS TO COMPUTE WATER SURFACE PROFILES 


pe M. Lara, A. M. ASCE and Kenneth B. Schroder,2 A. M. ASCE 


SYNOPSIS 


vo methods derived by the authors in their work with the Bureau of 
mation are presented for the computation of water surface profiles used 
ie development of tail water rating curves and backwater profiles. Ex- 
2s of each method are included which show the computational procedure. 


INTRODUCTION 


e derivation of water surface profiles serves two major purposes in the 
iu of Reclamation; one—to determine tail water rating curves, and two— 
ce the backwater curve above adam, Tail water curves are used in the 
n of power-plants, pumping plants, and energy dissipators such as still- 
sins. In the design of large dams, these curves also furnish useful 
nation for making stability and stress analyses. A primary use of the 
rater curve above a dam is for reservoir land acquisitions and 

1ents, Backwater information is also necessary in the design of new fa- 
S such as pumping plants, bridges, powerplants, or other structures 

d along or above the reservoir. 

e methods for computing water surface profiles developed herein are 
red to as Method A and Method B. Method A can be used in computing 
es for channels in which the flow distance between sections is equal for 
zments of the flow. Method B differs from Method A in that it is es- 

ly adaptable to conditions where the flow travel paths are longer for the 
channel than for the overbank portion. Both methods, however, are based 
roper hydraulic analyses of the conveyance capacities of each cross- 
nal subdivision. Eddy loss corrections and velocity head changes can 

e taken into account by either method. 


Discussion open until September 1, 1959. To extend the closing date one month, 
ritten request must be filed with the Executive Secretary, ASCE. Paper 1997 is 
t of the copyrighted Journal of the Hydraulics Division, Proceedings of the 
erican Society of Civil Engineers, Vol. 85, No. HY 4, April, 1959. 

ir. Engr., Bureau of Reclamation, Denver, Colo. 

ir. Engr., Bureau of Reclamation, Denver, Colo. 
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Theory 


The underlying theory of each method is based on Bernoulli’s Energy 
Equation which is repeated for convenience in the following discussion of it 
application to water surface profile computations. 

The conditions of an open channel reach in steady, non-uniform flow are 
diagrammatically illustrated in Fig. 1. The energy balance equation is wri 
ten: 


2, + a, + fe = 2 + qa, + an + hp + other losses 


Zo = streambed elevation referenced to a given datum (upstream sectior 


dg = depth of flow at upstream section 


hy, = velocity head of upstream section 
Z 1 = streambed elevation referenced to a given datum (downstream 
section) 


d; = depth of flow at downstream section 
le velocity head of downstream section 
hg = friction head loss 
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ENERGY OF OPEN-CHANNEL FLOW 


FIGURE 1 
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er losses = eddy, bend, and bridge pier losses. 


e two methods use basically the following hydraulic relationships to 
ate the various components of the energy equation. 


ye 


Q= AV, Vie V/A, by = 5 


a 


tt Ron) 
1.486 aR /? 


oe 


= discharge in cubic feet per second 

= cross-sectional area in square feet 

= cross-sectional velocity in feet per second 

= gravitational acceleration in feet per second squared 
= reach length between sections in feet 

= roughness coefficient representing the reach 


= hydraulic radius equal to cross-sectional area divided by the wetted 
perimeter 


is noted from above that Manning’s equation is used to compute hg. 
rally, Manning’s discharge equation is written as follows: 


2 1 
Q = Be AR /s Ss /e 
e S is the slope of the energy gradient. 


es 


1e total losses in the energy balance equation consist of the friction head, 
ldy loss, bend loss, and bridge pier loss. 

1e friction head can be determined by applying Manning’s equation in 

s of the friction slope, Sr: 


2 
Sete kina ae 
otal friction head for a reach then can be resolved by the following 
ion: 
Ris +S fe 
Hm ts ( 2 ) 


> Sf, and Sfo are the friction slopes computed by the preceding equation 
the representative hydraulic properties: n, A, and R for Sections 1 and 


simplified and practical concept of eddy loss can best be understood by 
ing the energy losses incident to the conversion of kinetic head to static 
and vice versa in open channel flow. While there are, at present, no 
ds for measuring these losses, it is known they do occur and manifest 
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themselves usually in the form of turbulence. Observations of channel flow 
where there are abrupt changes in velocity clearly display turbulence by the 
presence of boils and vortexes. This is particularly evident when there is | 
significant change in cross-sectional area. 

It has been a general practice in the Bureau of Reclamation to compute th 
eddy loss as a 10-per cent correction of the difference in velocity heads (Ah, 
when the static head is converted to kinetic head and 50 per cent for a rever 
condition. Thus, when the value of hy, - hyo is positive, a 10-per cent cor= 
rection is applied, and for negative values a 50-per cent correction is used, 

Departures from the 50-per cent value quoted above may be necessary, ; 
particularly, in cases where extreme changes in velocity occur in the channt 
reach. This condition results in transitions from a very narrow section to4 
wide section. Large standing waves are formed, which indicates the flow is 
bordering critical stage. As the flow progresses downstream, the energy 4 
dissipated in turbulence and upon reaching the downstream section, it reve 
to a more tranquil state. In some of these cases, the eddy loss could reason 
ably be increased to as much as 100 per cent of the change in velocity he 
By always applying a 50-per cent correction in the computations, a situati 
can result showing a higher water surface elevation at the downstream se 
(Z1 + dy > Zg + dg, an adverse water surface slope) which is seldom en- — 
countered in natural open channel flow. To correct this condition, an ad- ~ 
justment can be made by increasing the percentage applied. Eddy losses, at 
best, are only estimates based on the judgment of the hydraulic engineer. _ 

A search of the literature reveals that there have been no formulas de- 
veloped yet for a thorough determination of bend losses. A number of formi 
las resulting from various studies are listed on following page. : 


channel width 


: = hy from Yarnell and Woodward(1) 
inner radius 


A. hp = 0.21 


B 
C. hp = 0.38 hy by C. H. Yen(3) 
D 


0.0256 v2 x width of stream 


= i 4 
b graphic radius of curvature hy Pere ats Urick(4) 


“n” increased 0.001 for each 20° of curvature suggested by Scobey(2) 
3 

There are no set criteria which have been established to make.a practical ; 
evaluation of bend losses. Two methods are suggested for the determinatioi 
of these losses; first, they could be computed from one of the formulas lis b¢ 
above; and second, “n” values could be modified to include bend loss effects 
Ordinarily, such losses determined by the first method are very theoretical 
and require analyses usually beyond the limits of practical dimensions of th 
collected field data. It is recommended, therefore, that “n” coefficients be 
assigned values to include any significant effects of bend losses. | 
The hydraulic conditions at bridge crossings require careful analysis to 
determine bridge pier losses. King’s Handbook(5) recommends contracted 
and enlargement coefficients of 0.5 and 1.0, respectively, applied to the cha 
es in velocity heads at the approach, bridge, and emergent section. Nagler 
Yarnell,(7,8) and Johnson(9) have developed formulas for computing the hez 
loss through bridges which may be applied under certain conditions. 
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Methods 


| 
) 
: 
/ 
a 
is is a trial and error procedure which involves step-by-step compu- 
s. The method has a very wide scope of application and can be used in 
heral case involving steady, non-uniform flow. It is particularly appli- 
/ to the irregular channel in which the cross section consists of a main 
nel and separate overbank areas having individual “n” coefficients. Ve- 
7 head changes are taken into account by a weighting process and cor- 
ons can be included for eddy losses within the reach. Reach lengths 
2senting the flow path between sections, however, are assumed equal for 
nain channel and overbank areas. sg 
ferring to the diagrammatic sketch in Fig. 1, the objective is to de- 
ine the water surface elevation for a given discharge at Section 2 when it 
own at Section 1. 
hen the discharge exceeds main channel capacity, the flow also inundates 
.ood plain or overbank area. If channel bed conditions are different from 
» of the overbank areas, separate “n” values are assigned, and the cross 
on is subdivided for proper hydraulic analysis. 
ne areas and hydraulic radii are computed separately for each subdivision 
abulated as shown for Section 1 in Table 1. The conveyance capacity, Kg, 
| to 1.486/n AR2/3, also shown in this table, is computed for the main 
1el and overbank using the corresponding “n” value. Conveyance curves, 
3. elevation, are plotted for each section; an example is shown in Fig. 2. 
this method, the friction slope, Ss, is determined by the following 
‘Sis: 


2 Fl 
q = 286 pgi/ g'/2 (Manning) 
1.486 ..?/s 
Ky ae AR 
1/5 
Go Kae 
ng for Sr 
2 


Sola) 


1e Kg value used in the above equation is the total conveyance capacities 
2 main channel and overbank areas. 
.e discharges occurring in each subdivision of the cross section are de- 


ined as follows: 


/2 


Be 
Qin Bate 


: 
= total discharge in cfs 
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HYDRAULIC CHARACTERISTICS 


Saction No. l Main Channel 


n # 0.030 1.486/n = 49.63 


f 


Elevation A R R2/3 
5701 2 0.48 0.6 
5702 52 1.27 1.173 
57hls 4g 7.21 3.732 
5713 547 8.17 4.057 
* 5713.5 570 8.55 4.181 ; 
S.7r 4.253 


Overbank Area Ay 
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o ~ een x yo = a eV woe . eo - mn anew 


1000 10,000 100,000 1,000,000 
FIGURE 2 

Ss = friction slope 

Kg =total conveyance capacity 


) and Kg _ = corresponding elements for partial discharge and con- 
veyance capacity of the subdivision under consideration. 


ling a unit slope in the above equations, they are written thus: 
Q=K,x1 and Spe hal © 
ig one equation by the other 


3 for Qh 
QxK 


ae 


Q = 
a the preceding analysis, it is fies that the Kg value, in addition to its 
ion as conveyance capacity, can also be designated as the discharge oc- 
3 with a slope equal to unity. 
friction head, hg, is determined by averaging the slopes at Sections 1 
hen multiplying by the Aone as siatrated below: 


) 
| j Z { 
' | 
, > : 
: 
q 
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Se +5, 


5 2 
he = L( - ) 
The velocity head, h,, is derived by a weighting process using the partial 
discharges occurring in each subdivision of the cross section. This method 
is identical to the one used by the Corps of Engineers. (10) Velocities in eae) 


segment are computed by the basic equation, Vp = Qp/Ap> where Ap is the Y 
area of the particular segment. The weighted velocity head results from the 


following relationship: 
1 Q,) ' 


= oP) 


h 
V 


Eddy losses are computed from the formula M(hy; = hyg) where M is 0.5 
for negative values of hy, - hys and 0.10 for positive values. 
The following example iustestes an application of Method A using the 
of a fictitious stream, the Red Fox River, Colorado, for a discharge of 11,1 
cfs. The column numbers refer to the columns in Table 2. 


Step 1 


Enter the cross-section designation under Column 1. Column 2 shows # 
assumed water surface elevation. For Section 1, this has been established 
as 5714.0 from a predeveloped rating curve. The main channel and overba 
areas and conveyance values are entered in Columns 3 and 5. 


Step 2 


The friction slope, S¢,, is computed from the equation S¢ = (Q/Kg)2 as — 
shown in Column 6. Thus, in the example, this is equal to (11,100/129,1 
0.00755. Note that the total Kg (Column 5) is used. 


Step 3 


Column 9 shows the subdivisional discharges as computed by the equati 


Q, = 2Edp 
ae 


Substituting the values of the main channel in the example: 
s (11,100) (125,000) 

* TEER ea ea . 

A similar computation is made for the overbank area using its Ky value 0 

4,150. 


Step 4 


The velocities in Column 10 are obtained by dividing Column 9 by Colu 


for each subdivision. Velocity in the main channel equals 10,700 divided | 
594, : 


Step 5 


Column 11 is computed by squaring Column 10 and multiplying by Col i 
The values of each subdivision are totaled and entered. 


ve a1avl 
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Step 6 


The weighted velocity head, hy, is computed from the following relation- 
ship: 


Substituting in the example: 


which is entered in Column 12. 
Step 7 
The distance between sections is entered in Column 4. 


Step 8 


An elevation is assumed and entered in Column 2 for Section 2. Similar 
entries are made in Columns 3 and 5 as was done for Section 1 described in 
Step 1 above. 


Step 9 


The friction slope is computed as described under Step 2 and entered in 
Column 6. Column 7 shows the average of the consecutive slopes as follows 


0.00 A 

2200755 + 0:000°T5 . 9 ookor 
Column 8, the friction head, equals Column 7 multiplied by Column 4. 
Step 10 


Columns 9, 10, 11, and 12 are computed as described in Steps 3, 4, 5, an 
6. 


Step 11 


The algebraic difference in velocity heads is entered in Column 13. Intl 
example, this equals 4.85 - 0.43 = 4.42. 


Step 12 


The eddy loss Column 14 is taken as 10 per cent of Column 13 for positit 
values or 50 per cent for negative values. This equals (0.10)(4.42) = 0.44 ir 
the example. 


Step 13 


Column 15 is the summation of Columns 8 and 14. Column 16 equals 
Columns 15 plus 13. 


Step 14 


Column 17 is obtained by adding the water surface of the downstream _ 
section (Column 17—Section 1) to the value shown in Column 16. For exam 
this equals 5714.0 + 6.86 = 5720.86. When the result in Column 17 is equal 
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nn 2 to the nearest tenth foot, the computations are begun at the next 

Lg section. 

pd B 

lis method is also a trial and error procedure involving step compu- 

1s. It is used in cases involving steady, nonuniform flow. However, it 

#s from Method A in that the reach lengths representing the flow path be- 
sections are different for the main channel and overbank. The over- 
each length could be considerably shorter. Consideration is also given 

- variability of the hydraulic elements “n”, area, and hydraulic radius as 

jone in Method A. 

rain, the sketch in Fig. 1 representing the energy equation is basic to 

method. Knowing the elevation at Section 1, the water surface elevation 

ction 2 is assumed and the energy equation checked. 

ue determination of the friction head, hg, will be discussed first. This is 

nplished by applying Manning’s formula as follows: 


2 1 
q = 188 mp (9g /2 


sssing Sf as Z where L is the length between sections. Substituting 
*/e 
_ 1.486 ,,?/s ( Be ) 
cane L 


ituting Kg in the equation 
‘ing and solving for hg 


Li2 


yer, since there are two or more values for Kq/L1/ 2 computed because 
differences in flow paths and conveyance capacities for the main 
el and overbank areas, the summation () of Kg/L1/2 is used: 


Q 
Bian 
% ( L/2 ) 
he, + h¢ 
an friction head results from the equation h¢ = er er where hg; and 


re the friction heads computed for Sections 1 and 2, respectively. 
e discharges occurring in each sectional subdivision are computed from 
rmula listed below: 


a7 ( xf. / a n, /2 


locity heads and eddy losses are treated the same as in Method A. 
application of Method B is shown below using the data of a fictitious 
m, Silver Fox River for a discharge of 50,000 cfs. Column numbers 
to the columns in Table 3. — 
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Step 1 

Enter the pertinent data in the heading of Table 3. In the example, these 
are the Silver Fox River water surface profile computations for a discharg 
of 50,000 cfs. 
Step 2 


Enter the station number 209+40 under Column 1. 


Step 3 


Under Column 2, the main channel and overbank reach lengths between 
sections are entered. These are noted as 2,260 and 1,955 feet between 
Sections 209+40 and 186+80 for the main channel and overbank reaches. 


Step 4 


The known water surface elevation of 94.67 is entered in Column 3 for 
Section 209+40. Areas of each subdivision are entered in Column 4. 


Step 5 


Kg values of 973,000 and 604,000 are entered in Column 5 for the main” 
channel and overbank area, respectively. 


Step 6 


Column 6 is determined by dividing the individual values of Column 5 by 
the square root of Column 2. For example, the main channel value is com 


ed as follows: 
K 
a 000 
= ee = 20,500 
L Te (2,260) E ; 


’ 
Step 7 


The friction head, hg, in Column 7 is found by squaring the quantity of tl 
total discharge divided by the summation of the values in Column 6, as _ 


follows: 
Ki 00. 
a (e) 
ne «(9 safe) = Giiaoo) = 2 


Po ee see ee SS ee 


Step 8 


The subdivisional flows listed in Column 8 are computed by multiplyi : 
the individual values of Column 6 by the square root of Column 7. For ex 


anes the overbank discharge is computed as follows: 13,700 x 2.131/2 = 
ey ° ; 


Step 9 


Columns 9, 10, and 11 are determined as shown in Steps 4, 5, and 6 wid 
Method A. 


aq 
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he remainder of the columns are omitted since Section 209+40 is the be- 
mg section. The next step, therefore, consists of assuming the water 

ice elevation for the next upstream Section 186+80. This has been as- 

rd as 96.90 in the example. Similar computations are made through 

mn 11 as described for Section 209+40. The velocity head, hg, for Section 
80 computes to 0.45. The algebraic difference in velocity heads is 

‘ed in Column 12 as in the example hy, - hys = 0.75 - 0.45 = +0.30. 

: 


111 
dy losses in Column 13 are determined as in Step 12 of Method A. 


12 


olumn 14 is the mean friction head loss, hg, computed by averaging the 
2s in Column 9 as shown in the following computations: 


be 2.15 2 1.62 hoe. 


13 


olumn 15 is equal to the total of Columns 13 and 14. Column 16 is the 
yraic sum of Columns 12 and 15. 


14 


he water surface elevation in Column 17 is calculated by adding Column 
the preceding water surface elevation (Section 209+40) of Column 17. In 
xample, this amounts to 94.67 + 2.21 = 96.88. When a balance to the 

est tenth foot is obtained between Columns 17 and 3, the computations 
-omplete. 

ne next step involves a computation of the friction head, hg, between 

ons 186+80 and 148+10 (next upstream section) using the lengths of 3,870 
},265 as shown for the main channel and overbank. The water surface 
tion at Section 148+10 is then assumed and a new cycle of computations 
gun using the same process outlined above. 
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fRANSITION FROM LAMINAR TO TURBULENT FLOW IN A PIPE@ 


Closure by M. R. Carstens 


. R. CARSTENS, | A. M. ASCE.—The discussion(1) by Prof. Robertson is 
‘rough review of the paper in which valid questions were posed concerning 
nterpretation of the experimental results. 

he most important question of the discussion concerned the propriety of 

¢ the Blasius boundary layer in the discussion of stability of the temporal- 
veloping boundary layer ina pipe. Specifically, the question was raised 

| whether or not boundary-layer stability studies by Tatsumi(2) would be 

2 applicable. Neither the Blasius boundary layer nor the Tatsumi bounda- 
yer is the time developing boundary layer in pipe, since the Blasius 

dary layer applies to steady-flow over a flat plate without a pressure 

ient and the Tatsumi boundary layer applies to steady flow in a pipe in 
nlet region. In spite of the differences of origin, the Tatsumi boundary 

- and the Blasius boundary layer are remarkably similar, as shown in 
Iin which y is the distance from the boundary. Since the second deriva- 
of v with respect to y is important in stability analyses, velocity profiles 
appear to be quite similar as in Fig. I but result in radically different 

cal stability curves. Thus, further evidence of the remarkable similarity 
een the Tatsumi and Blasius boundary layers is shown in Fig. II in which 
itral stability curve for the Tatsumi boundary layer is superposed upon 
eutral stability curve of the Blasius boundary layer. Therefore, there 
ars to be no intrinsic merit in using the Tatsumi boundary layer for 
tative discussion of a case for which the stability analysis is not availa- 
On the other hand, the advantage of using the Blasius boundary layer for 
larison is that the stability analyses are much more comprehensive than 
ny other boundary-layer velocity profile. 

ne comparison of the stability of the temporally developing boundary layer 
ey can also be made with the stability of the Polhausen boundary 

(3) In Fig. III the plotted points were obtaineu from an analog computer 
ion of the Navier-Stokes equation with boundary conditions and physical 
ants of Run 29. The analog computer solution for the velocity profiles 

n in Fig. III is in the range of boundary-layer thickness for which tran- 

1 to turbulence was observed. The ratio of the displacement thickness, 

) the momentum thickness, 9, is about 2.43 for the four velocity profiles 
nin Fig. II. This value of 6*/O corresponds to a Polhausen shape 

neter of 3. Stability analysis of this Polhausen boundary layer indicates 
he critical or minimum value of the boundary-layer Reynolds number, 


roc. Paper 1450, December, 1957, by M. R. Carstens. 
rof. of Civ. Eng., Georgia Inst. of Technology, Atlanta, Ga. 
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Tatsumi (8"/r, = 0.0513) 
(8"/ry = 0.1056) 


Blasius 
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Blosius Boundary Layer 
Run 29. 3”/g = 0.034 
S/t = 0.047 
8" /r, = 0.056 
6"/rg =0.063 
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Figure III. Boundary-Layer Velocity Profiles of Run 29 


IRs, is 3500. (4) However, in all of the experimental runs, the inception fe) 
turbulence was always observed (see Table III of the paper 1450) at a valu 
of IR5» less than 3500. This observation tends to substantiate the discuss 
statement “there is some evidence that a timewise developing flow may be 
less stable than a steady flow.” 

In spite of the hazards of using a borrowed solution, that is, the stabilil 
of the Blasius boundary layer, the general nature of the boundary layer in: 
bility in the time developing boundary layer of the experiments is undoubt 
identical to that of the borrowed solution. The fact that the borrowed stab 
analysis offered a rational explanation of the observed transition to turbul 
is indicative to the author that the neutral curve of stability of the time de 
veloping boundary layer must be very similar, both quantitatively and qua 
tatively, to the neutral stability curve for the Blasius boundary layer. 

The discusser has correctly deduced the manner of obtaining the “poss 
neutral stability curves” which are shown in Fig. 10,.that is, the disturbai 
wave length was chosen which would pass the “possible neutral stability — 
curve” through the experimental point having the highest transition Reyno 
number. The hypothesis that the disturbance wave length, A, is proportio 
to the roughness height was based upon an analogy with the Karman vorte: 
trail behind a circular cylinder. In a stable Karman vortex trail, the spa 
between successive vortices is proportional to cylinder diameter and inde 
pendent of the flow variables. The analogy is suggestive that the disturba 
wave length, A, is proportional only to the roughness height or, in this ca 
the step height of the mismatched pipe boundaries at the junction. The m 
matching of the pipe ends could not be expected to be axisymmetric. Thu 
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Histurbance wave lengths, ranging downward in magnitude from the 
hum wave length at the point of greatest mismatching, could be expected 
finate at the pipe junction. Only the possible neutral stability curve for 
mgest disturbance wave length was plotted on Fig. 10. The other shorter 
jbance wave lengths originating at the pipe junction would have neutral 
jty curves lying below the plotted curve. The curve defining the limits 
tral stability of all disturbances would be bounded by a vertical line on 
ft in Fig. 10 at a value of IRs, of 420 with an upper boundary as shown in 
(0. The above explanation was the basis for the statement “With this 
pt there is no conflict in that some of the points of turbulence inception 
ated by the circular symbols on Fig. 10) are below the lower branch of 
irve shown.” 
= question was raised in the discussion concerning the amplification of 
urbance which is required to produce turbulence inception. The stability 
ses of small disturbances are merely indications as to whether the dis- 
ace will amplify or decay. The stability analyses do not indicate whether 
aplification of a disturbance is sufficient to result in transition to turbu- 
The presumption is that, when a disturbance reaches a certain magni- 
transition occurs. The order of magnitude of amplification required for 
tion is not a constant, as implied by the discusser, but must certainly 
d upon the initial magnitude of the disturbance. For example, an infini- 
al disturbance would require infinite amplification and a finite dis- 
ice, which is equal or greater than the limiting magnitude, would cause 
tion without amplification. Since the transition to turbulence was always 
ved some distance downstream from the pipe junction, the conclusion is 
1e disturbances were small (but certainly finite) disturbances that re- 
1 some amplification before transition occurred. 
e discussion concerning Fig. 7 appears to result from a misunderstand- 
\ brief review of the total transition occurrence will clarify the meaning 
;figure. Prior to opening of the downstream valve, the fluid is at rest 
thout the pipe. At the instant of valve opening a piezometric-pressure 
nt is applied to the fluid in the pipe. The unbalanced force system re- 
in acceleration. The developing flow pattern can be visualized as a 
1 core with uniform velocity distribution and a boundary layer in which 
locity varies from zero at the boundary to the core velocity. The core 
ty increases linearly with time during the initial stages of the motion. 
jundary layer is of zero thickness at t = 0 and increases with time. 
ow is laminar throughout. Meanwhile, small disturbances are being fed 
e boundary layer downstream from the pipe junctions. There is only a 
d range of core velocity, ve, and boundary-layer thickness, 6*, during 
these disturbances will be amplified as shown Figs. 9 and 10. During 
riod of amplification, a disturbance from the pipe junction is amplified 
turbulent burst is formed. Once a turbulence burst exists, large ad- 
mts in the laminar velocity profile are effected. Since the pipe flow is 
J in extent, these adjustments within the boundary layer also result in 
‘iable velocity changes in the central core adjacent to the turbulent 
Since the concept of a laminar central core riding over a turbulent 
iy layer appears to be untenable, the turbulence must penetrate into 
ntral core in some manner. Whether the velocity adjustment leads to 
lity in the central core or whether the disturbance of the turbulent 
is sufficient to cause instability and a radial spread of turbulence is not 
to the author. In any event the schematic drawing in Fig. 7, which is 
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based upon a study of the film strip data and pressure data as the downstré 
face of the turbulent spot passed through the pipe outlet, appears to be a 
rational picture of the downstream face of a turbulent spot which moves 
toward the pipe outlet. Essentially, the same type of downstream face ofa ~ 
turbulent spot is reported by Schubauer and Klebanoff 5) for flow over a flat. 
plate as follows: E 


“The turbulence is, in fact, transported downstream with the free- 
stream velocity and the lag at the surface is due to the time required 
for propagation inward through the slower moving air of the laminar 
layer. Probably the chief significance of the slower progress near “= 
the surface is that it gives rise to an overhanging leading edge.” , 
7 


Two errors in the paper were pointed out in the discussion. First, the 
number of wave lengths between the point of disturbance origin and the po 
of turbulence inception would be about 4 to 37 for turbulent spot 2 and abo 
13 to 31 for turbulent spot 3 based on wave lengths of 2.1 rp and 1.5 ro, re 
spectively. Second, the X/D value for spot 1 of Run 20 should be 12 rather 
than 112 in Table II]. A question mark in parenthesis should be placed in 
Table III adjacent to X/D value of spot 1 of Run 21. 

The author thanks Professor Robertson for his carefully prepared dis- 
cussion of the paper and thanks Messrs. Johnson and Meeks of the Georgia 
Tech Analog Computer Section for their help in obtaining the solutions sh 
in Fig. II. 
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| SOME EXPERIMENTS WITH EMERGENCY SIPHON SPILLWAYS2 


Discussion by Fred W. Blaisdell and Harold W. Humphreys 


RED W. BLAISDELL,! M. ASCE and HAROLD W. HUMPHREYS,2 A. M. 
=.—India, in addition to the countries listed by the author, has used large 
bn spillways for the close control of reservoir levels. But only India, to 
vriters’ knowledge, has used both “saddle” and “volute” siphons, 1) as 
‘call them. This discussion will be confined largely to the volute siphon 
various modifications of it. 

he volute siphon might be defined as a morning glory spillway which has 
sntrance covered by a dome. The dome is of greater diameter than the 
ance and extends below the crest elevation in the same manner as does 
op or hood of the saddle siphon. Also, like the saddle siphon, the volute 
jn uses siphon breakers which run to the top of the dome. In contrast to 
saddle siphon as detailed by Mr. McBirney, the volute siphon does not 

e the lower leg upstream or necessarily use a device to throw the nappe 
e€ crown, nor is it necessary to submerge the exit to insure priming. 
dian literature credits Shri V. Ganesh lyar with the discovery of the 

€ siphon in 1933 and its experimental development at Mandya in 1935-36 
y a Siphon 2.5 feet in diameter and 11 feet high. (1,2) Two volute siphons 
feet in diameter with an operating head of 43 feet were constructed at 
-onahally in 1938 as an experimental measure. Eleven siphons 18 feet in 
eter with a drop from the design head-water level to the invert at the 
way exit of 68 feet were built in 1947 at Hirebhasgar. 

is not necessary to use the dome to insure that the spillway will run full 
r suction. Examples of this may be found in the Heart Butte, North 

ta, and Shade Hill, South Dakota, spillways of the U. S. Bureau of Recla- 
on.(3) Ifa siphon is defined as a spillway which requires that the liquid 
ited above the headpool level, then the Dakota spillways do not qualify as 
ns. They do qualify as siphons if the definition is: 


closed conduit, a portion of which lies above the hydraulic grade line. 
iis results in a pressure less than atmospheric in that portion, and 
mce requires that a vacuum be created to start flow. A siphon util- 
es atmospheric pressure to effect or increase the flow of water 
rough it. 


roc. Paper 1807, October, 1958, by Warren B. McBirney. 
roject Supervisor, Agri. Research Service, U. S. Dept. of Agri., St. 
athony Falls Hydr. Lab., Minneapolis, Minn. 

ead, Hydr. Research Section, State Water Survey, Urbana, Ill. 
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The Dakota spillways qualify as siphons because, under full flow, the vertica 
legs are under pressures less than atmospheric and lie above the hydraulic 
grade line. 

The possibility that the pressures within these spillways may be so low a 
to cause cavitation should be recognized. The Hirebhasgar volute siphons 1 
in India are examples. These siphons consist of a dome-covered morning 
glory entrance 30 feet 8 inches in diameter which converges to a drop inlet 
shaft 18 feet in diameter, a 90 degree elbow, and a horizontal leg 18 feet in 
diameter. The drop in level between the crest and the exit center line is ab 
57 feet—considerably in excess of absolute zero pressure head even conside 
ing hydraulic losses. Cavitation damage should have been expected and did 
occur.\1) Cavitation damage also has occurred in the Marconahally siphons, 
This should have been anticipated since one-sixteenth size model tests indi- 
cate prototype pressures below absolute zero over much of the drop inlet an 
inside of the elbow. (5) 

In the Dakota spillways of the U. S. Bureau of Reclamation, the possibilit 
of extremely low pressures was recognized and means were taken to preven 
their occurrence. In the case of the Heart Butte spillway, a lip is used to ~ 
cause the stream to break away from the wall of the drop inlet. This insure 
atmospheric pressure at that point and that negative pressures in the drop — 
inlet will not exceed a predetermined magnitude. The drop through the Shac 
Hill spillway does not appear to be sufficient to cause cavitation problems. 

Full flow can be achieved for spillways of the morning glory type without 
special priming devices if the inlet is properly formed. The writers are of 
the opinion that the Hungry Horse spillway(3) would flow full if the head ove 
the crest is raised high enough. If this were to occur, serious cavitation — 
would be a likelihood, in view of the 487 foot drop through the spillway, unle 
means were taken to prevent it. Presumably the maximum head over the 
crest is sufficiently low and the barrel sufficiently large to insure part full 
flow at all times. It is important that designers recognize that morning glo 
spillways can flow full, in order to know what to design against if they want 
part full flow, or to know what to design for if they want full flow. 

The Dakota spillways required no priming or depriming aids other than 
shaping the inlet properly. The head-discharge curves for the Heart Butte 
spillway presented by Peterka(6) show that the flow through the spillway is 
controlled by a weir at low heads. At high heads, the control is by pipe floy 
when the inlet is submerged and the drop inlet is full. At the higher weir — 
heads, the falling water must partially seal off the drop inlet and the increa 
head must create a demand for more water than can be supplied over the 
crest weir at that headpool level. The excess demand is supplied by air en 
trainment and air flow through the spillway. This is, the writers believe, t 
partialization for which Mr. McBirney is looking. Partialization is achieve 
here without special devices. : 

Only relatively large spillways have been discussed up to this point. La 
closed conduit spillways designed to flow full, like the Heart Butte and Shac 
Hill structures, are few in number. However, the U. S. Soil Conservation 
Service (SCS) began building this type of spillway in relatively smaller size 
in the 1930’s. Closed conduit self-priming spillways built annually under : 
Public Law No. 46, 74th Congress, entitled “An Act to Provide for the Pro- 
tection of Land Resources Against Soil Erosion, and for other Purposes;” é 
Public Law 566, 83rd Congress, known as the Watershed Protection and Fl 
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téntion Act, number in the thousands(7) and many of them operate at de- 
} capacity each year. The hydraulic design of these SCS spillways is 
tical to the design methods which should be used for the relatively larger 
ways of the morning glory type. The senior writer has outlined the theory 
tiously(8) and recently has authored a series of papers\9 describing both 
ratory and field tests on a number of different inlet configurations. It 
ibe of interest to mention in passing that the hydraulics involved is that 
nt in elementary courses but that the application requires a thorough 
rledge of the conditions under which each of the elementary hydraulic laws 
applied. 
great many of the SCS-designed closed conduit spillways consist of a 
ical drop inlet from the base of which a barrel carries the flow through 
lam. The slope of the barrel may be either steeper or flatter than that of 
aydraulic grade line. If the drop inlet is properly proportioned, the spill- 
‘will flow full whatever the barrel slope may be even if the outlet is not 
aerged. Pressures within the barrel will be less than atmospheric and 
barrel will lie above the hydraulic grade line if the barrel slope is steeper 
the hydraulic grade line. The spillway then acts as a siphon if the broad 
uition presented earlier is accepted. 

series of tests on drop inlet closed conduit spillways that are pertinent 
is discussion was initiated in 1957 at the Urbana Hydraulic Laboratory of 
llinois State Water Survey under the direction of the junior writer. This 
way consists of a circular drop inlet five barrel diameters (5D) high by 
barrel diameters (5D/3) in diameter. The barrel diameter D is three 
es. The 30 per cent slope on which the barrel has been placed initially is 
iderably steeper than the hydraulic grade line slope. Water enters the 
inlet from the entire periphery of the crest. The head-discharge curve 
his spillway has the form shown in Fig. A. Section I of the head-discharge 
e represents weir control with the weir at the crest of the drop inlet. 
ion II is also weir control but air is carried through the spillway with the 
r. Section III represents pipe control. The head-discharge curve and the 
have been fully described elsewhere. (10) 
he primary purpose of the Illinois State Water Survey study is to de- 
line the performance of a flat plate anti-vortex device supported above 
‘rest and to develop design rules for the form and dimensions of the plate. 
initial flat plate is two drop inlet (riser) diameters (2 Dp) in diameter 
three guide vanes, 120° apart, attached to the plate. The guide vanes ex- 
radially outward from the outside edge of the drop inlet. The solid curve 
ig. A is for the case where the anti-vortex plate is located 0.556DR or 
2r above the inlet crest, that is, within the range of pipe control repre- 
2d by Section III. When the plate was lowered into the weir control ranges 
I, it was found that the plate sealed off the inlet as soon as the water 
reached the bottom edge of the plate. This prevented access of air to 
lrop inlet and the flow caused the creation of a partial vacuum. The 
um sucked air under the plate in sufficient quantity to satisfy the demand. 
3asing the water flow caused little or no change in headpool level until 
pillway flowed full of water alone. Air was sucked into the spillway in 
ints to maintain a vacuum just sufficient to permit the water flow to take 
2 at constant headpool levels. The dash lines in Fig. A show the head- 
large curves for successively lower elevations of the anti-vortex plate. 
headpool level in essentially steady, exhibiting none of the fluctuations 
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tn by Mr. McBirney’s Fig. 4. Is this not the partialization and steady 
pool level that Mr. McBirney is looking for? The possibility of positively 
electing the priming head within certain presently undefined limits by 

ng the flat plate at the correct height above the crest should not be over- 
ed. 

0 give actual figures for comparison, Mr. McBirney’s spillway apparently 
a throat area of 1 x 0.5 = 0.5 square feet. The linois spillway has a 

el area of (0.25)2 7/4 = 0.0491 square feet. The linear multiplication 

or to secure comparable priming heads is then yV'9.5/0.0491 = 3.19. Plate 
nts as low as 5/8 inch caused the llinois spillway to prime. Increasing 
ito the scale of Mr. McBirney’s tests, the priming head becomes 0.17 feet. 
corresponds to Mr. McBirney’s full flow priming head of 0.52 feet for 
standard siphon design and 0.42 feet for the proposed design. 

jroper positioning of the flat plate anti-vortex plate used as a siphon 

1er results in a single-valued head-discharge curve. The curves of two 

Ss shown by Mr. McBirney in Figs. 6 and 11 for the standard and proposed 
zns do not exist for the flat plate priming device in the positions shown in 
A. However, two loops can be obtained for the drop inlet entrance by ex- 
ing skirts from the periphery of the cover plate down below the crest of 
weir. The spillway would then prime on rising stages at the same head as 
e skirts had not been used, but the prime would be held until the pool level 
hed the bottom of the skirts, or to the siphonbreaker if one is used. The 
ted arrangement is, of course, what is called a vortex siphon spillway in 
Le 

he priming time, as used by Mr. McBirney for the saddle siphons, does 
sxist for the drop inlet closed conduit spillway. Priming for the latter 
way takes place gradually as the headpool level rises and does not occur 
tively suddenly as for the saddle siphon. 

he gradual increase and decrease of flow obtained with the drop inlet 

on also reduces the erosion damage downstream from the outlet which 
McBirney says results at many locations from the prime-break-prime 

e of intermittent operation. 

he writers do not say that the drop inlet closed conduit spillway, either 
gut or with the flat plate anti-vortex device and primer, is a complete 
tion to all siphon problems. They do maintain, however, that the potenti- 
2s of the drop inlet siphon should not be overlooked. 
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SYNTHETIC FLOOD FREQUENCY@ 


Discussion by Max A. Kohler 


MAX A, KOHLER, |! M. ASCE.—Mr. Snyder’s approach to generalized, or 
thetic, flood frequencies can best be described as an elaboration and im- 
ivement of the so-called “rational method.” It involves the assumption that 
‘n-year flood can be derived from the n-year rainfall even though the two 
ints do not necessarily bear the relation of cause and effect. 

Specifically, the proposed method involves (1) determination of the n-year 
nfall for the basin with specified area and concentration time Tg, (2) con- 
sion of n-year rainfall to n-year runoff volume and (3) determination of 
flood peak corresponding to the runoff volume. This discussion concerns 
marily steps (1) and (2) of the procedure since it is believed that Snyder’s 
any of several other available techniques are adequate for converting run- 
volume to corresponding peak flow. 

Before taking up specific aspects of Snyder’s procedure, it is perhaps best 
‘onsider the basic underlying assumption that a simple relation exists be- 
en the n-year flood peak and the n-year rainfall. From a probability point 
riew, this implies that the composite probability for all factors which af- 

_ the relation between rainfall amount and peak flow is unity. Therefore, 
only direct means of relation n-year rainfall and flood peaks is through 
istical analysis of n-year data. This means that flood-frequency data 

st be used in deriving the generalized relations. There is no indication on 
part of the author that flood-frequency data are necessary or were used. 
tainly there was some basis for deciding that the required relation should 
> “slightly higher runoff than the average” (indicated by average monthly 
ifall and runoff). Wouldn’t one expect a different rainfall-runoff relation 
small basins (Tc < 4 hrs) where the major floods result from summer 
wderstorms than for the larger basins where such floods result from gener- 
ains on wet soil? It is believed the author has tended to account for this 
lation by an unrealistic reduction of point rainfall for T, and drainage 

2. Thus he chose to derive an area-reduction relation from storm data in 
ference to using one properly derived from frequency data.(1) It is not 

ir why his results indicate a flatter curve than one based on frequency 

1, while just the opposite would be expected. 

f it is visualized that Snyder’s procedure constitutes a method of generaliz- 
or regionalizing, on the basis of existing flood-frequency data, certain 
ifications would appear to be appropriate. Although not absolutely neces- 
y, it is suggested that the empirical tie between n-year rainfall and flood- 
< be confined to the rainfall-runoff relation. In this case, the n-year 


Proc. Paper 1808, October, 1958, by Franklin F. Snyder. 
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runoff can be computed from Snyder’s relations for all basins with establishe 
flood-frequency curves. These values can then be related to corresponding 
values of n-year basin rainfall derived from the point-rainfall frequency re= 
lation (author’s Fig. 2) and one of areal rainfall vs duration (T,) and point 
rainfall.(1) This is the manner in which the accompanying chart was de- 
veloped. 

The values of peak discharge (Qp) for various basins and return periods 
were obtained by Gumbel analysis with appropriate adjustment to give partia 
duration data. The average rate of direct runoff (I,) for time Tc was obtaine 
by use of the author’s formula, I, = Qp/ 500A, where A is the basin area in © 
Square miles. 

The first part of the accompanying relation, namely, the area curves in 
upper left, was derived from the 10-year rainfall curve of the author’s Fig. 
The 10-year rainfall amounts for various concentration times, or durations, 
were first plotted and the zero-area curve was drawn through the points. 
curves for other sizes of area were then constructed on the basis of the dept 
area relation(1) given in U. S. Weather Bureau Technical Paper 29. Althougi 
the depth-area relation is based on rather limited data, it is considered to 
the best available for the purpose. The second part of the writer’s relation, 
i.e., the return-period curves in the lower left, was obtained by entering the 
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adrant with T,, going horizontally to the zero-area curve, then down- 
o the appropriate rainfall intensity (indicated by the author’s Fig. 2). 
point so determined was labelled with the corresponding return period 
{smooth family of curves were fitted to the points, (2) 
should be pointed out that the return-period curves cannot be drawn to 
iy fit areas of all sizes, but the errors are insignificant with respect to 
rrors of the procedure as a whole. The final curve family (Te curves) 
serived by plotting n-year basin rainfall intensity computed from the first 
uadrants against corresponding values of n-year runoff intensity I, and 
mg the points with Tc. 
e accompanying table presents the basic data for eleven natural basins 
general vicinity of Washington, D. C., which were used in deriving the 
rves. These basins were selected because flood-frequency curves based 
records of observed and synthesized(3) annual peaks were available 
ire unit hydrographs for estimating Te. 
ie application of the writer’s relation for estimating n-year runoff in- 
ry is straightforward and needs no further amplification here. The pro- 
es described by Snyder for obtaining T, and for converting I, to the n- 
peak flow are equally applicable when using the writer’s accompanying 
on. 
conclusion, the writer wishes to commend the author for presenting new 
timulating ideas on an important phase of hydrologic design. 
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IDIVINING RODS VERSUS HYDROLOGIC DATA AND RESEARCH@ 
Discussions by Erhard E. Dittbrenner and V. M. Yevdjevich 


RHARD E. DITTBRENNER, | M. ASCE.—Mr. Langbein presents a scholar- 
i comprehensive review of the status of the most critical need in hy- 
Sy. It would be difficult to add anything useful to his presentation on hy- 
sy. However, when he includes “waterwitching” in hydrology—even if 
aS a strawman to be knocked down and used as a reference—it is time to 
rate fact from fancy. 
r. Langbein states: “Modern hydrology offers no support of water 
ery.” Ina certain limited sense this statement may be accurate enough. 
irther question: “. .. can we be scientific about witchery,” is obviously 
rical, as his implication is that since witchery is perhaps in the realm 
ery, it cannot be scientific. 
ie answer to his question is, however, an unqualified “Yes.” We can be 
tific about “witchery.” To do so, however, we must go to that portion of 
cal science which has to do with electrical phenomena. This might seem 
e it out of that portion of science which we have artifically labeled “hy- 
zy.” It is not important to argue the point, but it seems pertinent to 

out that as we push back the frontiers of our knowledge in our artificially 
ribed categories of science it becomes clear that Nature has no such 
"ary divisions. More and more evidence substantiates a belief perhaps 
ancient than water witching—that there is a oneness or unity in the 
rse, 
e writer happens to be a “water witch” or dowzer. Many years ago he 
nvited to try his hand by a young dowzer who suspected skepticism of the 
‘my attitude. The dowzer pointed out that many people have the dowzer’s 
ilities without realizing it. One trial demonstrated the truth of that 
nent and the basic reason why water witchery has persisted over the 

and through the ages. Once a person experiences the phenomena which 
the switch down, he will never be bothered by the scoffing of the “scien- 
mind. In the writer’s case, the switch broke in his hands on his first 

It has never done so since, but since he never worked at the art very 
this may not be very significant. 
ce there had to be some scientific explanation for the phenomena, search 
rature and query among those more qualified in science than himself 
7 brought some answers. The “dry” holes can be explained by the 
edge that the forces involved in the phenomena are obviously too small 
ak the twig, but act on certain muscles which are thereby impelled to 
about the action of the switch without volition on the part of the dowzer. 


oc. Paper 1809, October, 1958, by W. B. Langbein. 
dr. Engr., U. S. Bur. of Public Rds., Albany, N. Y. 
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The forces involved vary constantly—during the day, with the weather, and 
from season to season. Some experienced dowzers are aware of this. Oth 
obviously are not. Since the dowzers are not aware of the source of the pi 
nomena or their powers, they are unable to explain these variations. 

Two scientists in England have studied this phenomena extensively. Maby 
and Franklin, in “The Physics of the Divining Rod,” London, 1939, cover the 
subject comprehensively, but by no means exhaustively. The phenomena has 
been measured. Instruments have been devised to determine quantitatively 
the field in which the dowzer works to make his determinations. People hay 
been checked to determine their powers of “divination.” About one person i 
ten has these powers in some degree. Their investigation did not extend to 
the point of determining the reasons why some do and others do not have 
them. So far as this writer knows, no one has made any effort in this di- 
rection—in fact, since Maby and Franklin’s publication, no further scientific 
literature of that nature on the subject has appeared to the writer’s knowled 

In the writer’s own experiences, a routine laboratory experiment during 
college required the measurement of each student’s bodily electrical re- 
sistance. His resistance was way out of line with that of any other person if 
the group. The notes are gone but the memory says it was lower—perhaps 4 
third or a fourth of the “normal.” Some interim reading—also not then note 
—indicates that persons who become alcoholics must be more careful aroun 
electricity than others, as they are more easily electrocuted. There is no 
implication that persons with low electrical resistance tend toward or becom 
alcoholic victims, but there is an implication that bodily resistance to elect 
currents vary greatly and these variations have significant physical and 
physiological results. Medical research some 25 or 30 years ago, has es= 
tablished the fact that the human body is a generator as well as a conductor 
electricity. It is logical that generation and restivity would bring about a © 
difference in the electrical field about the human body and its reaction with 
the field about water “veins.” : 

Here, physiology and the physical sciences appear to have a common int 
est. So far, nothing has gone beyond what Maby and Franklin have reportet 
so far as the writer knows. It would seem to be an interesting field for bot 
the medical and the earth scientist to explore. 

The question now arises—“so what?” Granted that the dowzer has defin 
powers, is he useful to the geologist or the hydrologist? Mr. Langbein 
himself answers that question. In the absence of completely detailed geolo 
cal and hydrologic subsurface knowledge, the dowzer certainly can be usel 
In the writer’s opinion, the likelihood of any other knowledge supplanting hi 
appears remote. So far as the writer’s knowledge of geology and hydrolog 
permit, neither of the latter branches of science will produce the detailed 
knowledge available to the informed and experienced dowzer. It seems p : 
ble that an extension of the Maby and Franklin research could enhance the 
dowzer’s powers, and mayhap even elevate him to a “professional” catego! 

This would appear a subject of useful and profitable research. There 
certainly is no further basis for scoffing. 


V. M. YEVDJEVICH,! M. ASCE.—Two very interesting problems are d 
cussed by Mr. W. B. te a) contrast of the divining rods as an 


1. Hydr. Engr., Foreign senior scientist, National Bureau of Standards, 
Washington, D. C. 
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Ogic decision-maker versus the modern basic-data programs; b) con- 
of the increase of gathering basic-data versus the deficiencies in re- 
h necessary to guide, analyze and interpret the basic data, and to obtain 
ne general principles and the regional hydrologic characteristics. 
is clear that modern hydrology does not support “forked stick” pro- 
2 which does not have a serious scientific background. But the hydrolo- 
‘a be aware of some facts which contribute to the persistence of the 
*s activities, so that the hydrologist may overcome his seemingly ap- 
- disadvantages. Ground water hydrology, frequently very complex, im- 
/Some uncertainty in decisions upon the geohydrologic specialist, while 
wser’s short cut way to decisions makes a spectacular impression. The 
ir is often a man well acquainted with local conditions, and makes de- 
is by analogy, based on experiences. He is sometimes a kind of living 
sasy for a farmer to consult. He is available, and what is important, 
These are the reasons why the use of the divining rod is so persistent 
ill continue to be so for relatively smaller activities until replaced by 
sngineering procedures, competitive from other points of view. 
2 collection of basic data in hydrology has been started and made under 
general assumptions for the expected water resources development and 
1. The basic concepts laid for the coliection of hydrologic data depend 
he future water resources development. The first concepts were based 
expected individual and mostly single purpose projects. Next step in 
pment was the multiple purpose use of water resources. The future 
robably ask the total and very integrated water resources development 
tate of great demands for water control and use. The basic concepts 
ning the collecting of hydrologic data today should anticipate the situation 
er resources development 2 - 3 decades or more from now. This is ~ 
‘r reason to support Mr. Langbein’s emphasis for a better scientific 
ition in general collection and interpretation of basic-data. 
en a demand for higher scientific efforts in hydrology is forwarded, the 
sary problems connected with it have to be emphasized. The efficiency 
research in hydrology is an important question. Surveying the hydro- 
iterature of the countries with relatively abundant funds and large 
r of people employed in research, the impression is that the results do 
rrespond to the efforts. The quality of the results and the possibilities 
2 useful applications are rather limited in comparison with the quantity 
lished papers and other literature. The greatest results for given funds 
eir best application to the current needs and problems can be achieved 
y timely and thorough preparations for expending research. 
» timely solution of many problems can help the rapid development of 
“ch in hydrology, as for example: the systematic preparation of highly 
ed engineers and scientists in hydrology; the development of good pro- 
‘s in putting forward and in selecting the topics to be studied; the ex- 
: of experiences, ideas and new methods; the type of work, in teams, 
lual or both; etc. 
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QUEUING THEORY AND WATER STORAGE 


Discussion by B. W. Gould 


by GOULD.!—The author’s paper shows that the techniques of mathe- 
statistics may be applied in the solution of the problem of determining 
red water storage, and also in estimating the probable behaviour of a 
yoir already constructed. 

2 service function: 


D=b+kS 


e interpreted in several different ways, each of which gives a different 


avoid confusion, the writer has adopted the following notation: 


a= Average rate of discharge in the interval between time x and time 
y. 
hs Average rate of inflow in the interval between time x and time y. 
= Instantaneous rate of discharge at time x. 


= Volume of stored water at time x. 


> unit of storage used by the writer is the same as that used by the 

, being equal to unit flow for unit time. 

» first interpretation of the service function to be considered is that the 
rge during an interval of time is determined by the storage at the end 
interval. This type of control would be difficult to achieve in practice, 
considered because it gives the same solution as that obtained by the 


Dy xe1 9 = D+ Whee 
milarly 
Det es te koe 


sic equations obtained by the author may be verified by subtraction, and 
ution of 


D 


ext for (S44 -. S, 


I, 9x1 
ther interpretation of the service function could be that the instantane- . 
scharge is dependent on the instantaneous value of storage, i.e. 


1c. Paper 1811, October, 1958, by W. B. Langbein. 
r., Horsham Waterworks Trust, Vic., Aust., formerly lecturer in Civ. 
., Univ. of New South Wales, Aust. 
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and 


Deay FD + kB Ly y: 


x+1 
By subtraction and substitution, this gives 
D+ = dD, = k(1y x41 = Dy yx+1 ) 


If the time interval is sufficiently small, Dx, x41 is approximately equal 
1/2 (Dx + Dx41). Substituting this value and rearranging gives: 


2k raped Oe 
D = =a [| oS ———— D 
x+1 eo ae ae Xy x+1 T° ey x 
instead of the author’s 
Fe k 1 
2 atk. Beadete sie taeeed 


In yet another case, it may be taken that the discharge could be determi 
by the amount of water in storage at the beginning of an interval. This type 
of control would be suitable for irrigation of an area containing both semi- 
permanent (e.g. citrus or vines) and annual (e.g. market garden) developms 
where a certain minimum amount of water is required to maintain the sem 
permanent vegetation, whilst the area of annuals to be planted in any year 1 
be varied according to the amount of stored water available at the beginnin 
of the year. This gives the equations: 

D,. xt! b + kS, 
and 


Dy-1,x° = b + kSx_4 


Subtracting, and substituting for (Sy - Sx-1) gives ‘ 


Oz EE ee Deen 


An estimate of the discharge in the nth year, given by successive substi 
tution of “D” values, is: . 


D anne 
n=15n 2yn-1 * C-K)Tn 3nd + (1-k)?T, yn 3 


ete. + -K§t-1) 
(1-15 Cre iano Cotes 


From this, the standard deviation of the estimate of the discharge can be- 
calculated as. 


Asoo was pes ieee . 


instead of . . 
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wen by the author in Eq. (6). 
}comparison between the author’s Eqs. (12), (13), and (14), and the corre- 
Hing equations derived by the writer for the case where the discharge 


ig a period is based on the storage available at the beginning of the peri- 
5 given in the following table: 


nn. Discharge based on Discharge based on 
De final storage. initial storage. 
| aS 2 
Dd {ec 2 ba -b é k = — 
(Z-m)? - Cv)? (-m)* + (K-b) 
F ee ne = eb) 5 = Gem)? + (Kv)? 
(X - b) (KX - b) 

a a z & 
MY Set v Ta gag yt Ny, 

4, Td. 


is to be noted that where the discharge is determined by the amount of 
‘d water at the beginning of the period under consideration, the equations 
k = 1 whenb=m. In this case, any flow in excess of the minimum flow 
y year may be considered as being stored for use in the immediately 
wing year, thus giving time to make plans for its utilization. Under such 
tions, the storage required to provide 100 per cent utilization would be 
ifference between the maximum probable inflow, and the minimum proba- 
atural inflow. This is verified by the equation which gives the required 
ge as 2to. 

3 in the case considered by the author, the storage required increases 

ly as b approaches x, and becomes infinite when b = X. 

1e accuracy of the “probability routing” methods described by the author 
amples 2 and 3 in the paper, is limited only by the number of storage 
vals used, the number of successive approximations made, and the ever- 
nt errors of sampling. Electronic computation using standard 

“ammes for the solution of linear simultaneous equations has been used 

© writer, who found that smaller intervals of storage volume could be 
and that the need for successive approximations was eliminated because 
jlution was provided in a direct manner. 
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